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ABSTRACT 


Sideslip and roll responses of an aircraft are coupled 
calling upon the pilot to use considerahle skill to control th 
aircraft. The requirement that pilot's joh he made easier, 
enabling him, in addition, to achieve better flight-path 
accuracy, motivates the need to decouple these lateral outputs 

decoupling is achieved by diagonalising the system 
transfer function matrix of the lateral set of equations 
derived, using state feedback. The decoupling problem stated 
involves the problem of pole-assignment to achieve desired 
stability characteristic. 

The aircraft data made use of is that of an airplane 
similar to X-15, a high-performance aircraft. 

Decoupling controller involving state feedback has been 
designed for this aircraft. The relative stability of the 
lateral dynamics of the aircraft has been improved by pole- 
assignment in order to conform to military handling qualities 
specifications. 



CHAPTER I 


IHTRODUGTIOH 

In the crowded terminal areas of today's air transportation 
and for the militaiy requirement of formations of flight during 
a mission, it v/ould be helpful both to provide a more precise 
flight --path control and to reduce pilot workload. These two 
improvements which are often conflicting, can be achieved by 
providing a decoupled flight control system. The core of the 
decoupling concept is that for a multivariable system - the 
aircraft in this case »• each control should affect a single 
output without disturbing others. Pilot workload is reduced 
since coordination of controls is no longer required. Addi- 
tionally, flight-path precision can be enhanced since the pilot 
can focus his attention on a particular output. 

It can be appreciated that after decoupling has been 
achieved, it should be ensured that our system remains stable 
in its dynamics. Because of the current design trends tox^ards 
the use of low aspect ratio, sweepback and higher wing loading 
which form some of the characteristics of high-performance 
aircrafts, stability of lateral dynamics of the aircraft has 
bocome very important. This calls for consideration of 
motions outside the' plane of symmetry, the airplane being free 
to roll about its x-axis, yaw about its z-axis and translate 
along its y-axis. The pilot effects these motions by causing 
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angular deflections of control surfaces. Since the lateral 
outputs so obtained, broadly specified as the bank angle and 
the sideslip angle, are often interacting, the pilot is 
called upon to acquire considerable skill in order to simulta- 
neously manipulate the inputs and successfully control the 
aircraft. The problem of decoupling and stability must, 
therefore, be viewed together. 

The decoupling problem with its application to aircraft 
lateral control systems h8,s received some attention and 
several approaches are possible. One method due to Hall [ 5 , 9 ] 
and Yore [5,10] uses model following technique. They specified 
a model that was decoupled and had desirable dynamic characteri 
sties, and then chose a controller so that the original system 
would follow the model. Though it is a simple approach, it 
does not always work. In fact, the controller used by Hall 
minimizes a certain norm but need not make the norm zero. Thus 
perfect model following is not assured. Por Hall's aircraft 
if the tangent of a certain reference variable is not zero, the 
method fails. 

Another approach taken by Montgomary and Hatch [ 7 , 5 ], 
called differential synthesis, is without theory and allows 
direct calculations of the set of feedback gains and control 
interconnects that yield arbitrarily selected flying qualities 
parameters. They specified that the required elements in the 
output/input transfer function of the aircraft be zero and use 
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a clever computational method to solve the decoupling 
problem. 

In yet another approach. Cliff and lutze [5j applied the 
Geometric Decoupling Theory developed by Wonham and Morse, 
following this geometric approach, they first found the most 
general feedback and control interconnect matrices and then 
used a quadratic penalty function method to arrive at specified 
handling qualities. Though this was quite a flexible method, 
it presented a numerical difficulty in the Gram-Schmidt 
orthonormalisation procedure which it involved. This difficulty, 
which necessitated recognition of a certain vector becoming zero 
during computation, if not overcome, produces erroneous results. 

In this work we follow the approach of Palb and Wolovich 
[2] and Gilbert [l]. The idea is to find feedback and control 
interconnect matrices such that the output/ input transfer fun- 
ction is diagonal and nonsingular. This method, making use of 
state feedback, also provides a means of achieving a specified 
pole -placement while simultaneously decoupling the system thus 
resulting in accomplishment of a stable decoupled system. 

Apart from the relative ease of solution of pole-pla.cement, 
Gilbert’ s approach, should in principle form the basis for a 
useful synthesis procedure as his theoretical results' were 
proved by constructive arguments and provided complete structure 
of the decoupling solution. His methodology was reduced to an 
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algorithmic form [3,4] and then mechaniaed by a general 
purpose computer program which yields all the necessary 
structural results and synthesis data. 

Chapter II of this work derives, ab initio, the lateral 
directional eq\xations of motion of the aircraft in the 
state-vector form. Variables of interests and the inputs are 
defined therein. Chapter III states important results of the 
decoupling theory, develops the algorithm, explains the steps 
involved in the computer program and states the decoupling 
problem. Chapter IV deals with the practical application and 
Chapter V reports the conclusions stemming from the study. 
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CHAPTER II 

MTERAL jjYHAMTCF! OP AiaGRi.?T 

Tlie purpose of this chepter is to derive the equations of 
la.teral motion for a rigid airplane in the conventional form. 

In s.c doing, several a.ssumptions will have to he made;. Tnese 
assumptions and their implications will he pointed out. The 
conventions,! form of the equations will then he trsaisformed 
into state equations hy defining suitable state variables. By 
specifying the output variables of interest, the output vector 
equation will then he formed. A brief mention is made 8,bout 
measuring of state variables. 

2.1 GOORBIMTE SYSTEMS ARB EXTBRIAI PORCES 

Consider Pig. 2.1, whore the airplane is shown as a bodj'' 
in space. To keep track of its motions, an inertial coordinate 
system X’ Y’ V is introduced. In most of the problems, 
earth's motion is sufficiently slow to enable us to consider 
X' Y' Z' to be earth fixed [21.]. 

The aircraft is a continuum of mass particles dm with 
position vector r' in the X' Y' Z' coordinate system. Each 
mass particle is subject to the force of gravity. The 
force per imit volume is. 


( 2 . 1 ) 
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2.1 Coord inatejjfgt^ 
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where is the mass density of the airplane and g,, the 
acceleration due to gravity. Observe that g is always 
oriented along earth's Z' axis. 

Mass pa.rticles located at the surface of the aircraft are 
also subject to an aerodsmamic force F (including thrust) per 
unit area. These two external force systems (i.t, R and F) 
are assumed to be the only ones acting on the airplane. 

2.2 EQUATIONS OF MOTION 


An aircraft has six degrees of freedom : three transls-- 
tional and three rotational. Conservation of both linear and 
angular momentxm, in terms of Newton's Second law, yields 
the following equations of motion; 


di!* * 

It I \ at '^A 


/ S. dV + J F dS 
s 


( 2 . 2 ) 


(Rate of change of 
linear Momentum) 


(Applied forces) 


d 


dr' 


dt ^ ~ Pa dt 

(Rate of cha.nge of 
angular Momentum) 


dV =/ r ' X g dV +/ r' X F dS (2.3) 


V s 

(Applied moments) 


where / dV and J dS represent volume and surface integrals 
V s 

respectively . 


The total mass of the airplane is 


m = / P. dV 

V ^ 


(2.4) 
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We assiime that the aircraft mass m is constant, that is, 



This assiimption is justified as the change of mass is 
small over periods of time used in dynamic stabilil^r calcula- 
tions. 

It is also assumed that' the mass distribution is constant 
with time. Effects such as fuel slosh are neglected. 

Now a new coordinate system ^ T Z, having its origin at P, 
the centre of mass of the airplane is introduced. Point P is 
kept track of in inertial space by means of the position 
vector r^. All mass particles are referred to P in X Y Z by 
means of position vectors r. X Y Z is called a body-fixed 
coordinate system. It is rigidly attached to the aircraft and 
thus moves with it. However, its orientation rela^tive to the 
airplane remains free to be selected. Clearly, for any mass 
particle, 


£'= Ep + £ (2.6) 

Since P is the centre of mass of the aircraft, 

/ r P. dV = 0 


(2.7) 
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This leads to the followinssr definition for r' 


< = i ^ 37’ dV 

~P m A - 


( 2 . 8 ) 


The left hand side of eqn. (2.2) can, therefore, be rewritten 


P f. d7 = ^ ~r I P. (r'+ r)dY = fr m r' 

dt ^ A dt dt dt A p — ^ dc dt —p 


j ar ■ 
Cl *T) 


m TT VoT^ = 21 
dt dt 


(2,9) 


where 


is defined as the velocity of the 


airplane centre of mass 


The right hand side of eqn. (2.2) can be written as 


/ P a S. + / H dS = mg. + P. + 


■A ^ 


where is the total aerodynamic force and Pgi is the total 
thrust force. Sqn. (2.2) would, therefore, become 


“ d t = + 4 + ^T 


( 2 . 10 ) 


Substituting eqn. (2.6) into (2.5) and rearranging , 


d f 

dt ^ dt Pj. dV = / r X P dS 

V s 


( 2.11 ) 


Pqn. (2.11) can be written ai 


dt / - ^ dt “^A 4 ^ 


( 2 . 12 ) 
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where is the total aerodynamic moment vector and ^ is 
the total thrust moment vector; 

The volume integral on the left hand side of eqn. 

(2^12) is a time dependent function. The time dependence 
can be eliminated by a. switch in the coordinate system. By 
rewriting eqns; (2,10) and (2.12) with respect to system 
X Y Z instead of X' Y' Z' , it turns out that the time- 
dependence of integration is eliminated. However, system 
X Y Z is a rotating coordinate system. From classical 
mechanics [21 ], it is well Icnown that when a vector A is 
transformed from a fixed to a rotating coordinate system, 
the following relation holds, 

gri = 3t 0. X A (2.13) 

Fixed(X'Y'Z' ) Rotating (X Y Z) 

where w is the angular velocity of system X Y Z with respect 
to system X' Y* Z’ and is obviously identified as the angular 
velocity of the aircraft since system X Y Z is rigidly 
attached to it. I 

Applying eqn. (2.13) to left hand side of (2.10) gives, i 

I 

dV ^ 

dt ~ ^ + (2.14) I 

I" 

I 

[ 

1 

I 



Therefore, 


m(ip + w Ip) = 


(2.15) 


Similarly, the left hand side of eqn. (2.12) can he written as, 


dt 


i r 


X 


dr 

dt 




dV = / r X 

V 


d^ 

dt 


dt 


A 



X 


d_ 

dt 


(r + a> 


X 


= / E 


X 


(r+wxr+2wxp + aj 


Uow, as the aircraft is assumed to he rigid, 
eqn, (2.12) yields, 


dV 

(‘f L) 

( 2 . 16 ) 

r = r = 0. 30^ 


/ ^ ^ (a ^ L + <iL ^ (l^ ^ £)) dV = M. + Mm (2.17) 

Y —A —'i 

The only time derivative occuring in eqn. (2.17) is “ , 
the angular acceleration of the aircraft, which is independent 
of the volume integration. In other words, the time-dependent 
integral has indeed been eliminated. Eqns. (2.15) and (2.17) 
are the vector form equations of motion written in the hody- 
fixed I Y Z system. 


The equations of motion (eqns. (2.15) and (2.17)) are 
now to he expressed in terms of their components along the 
coordinate directions of the X Y z system. The coordinate 
directions are given hy the unit vectors i,j and h, in the 
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usual manner. See Table 2,1 in this connection. 

2.2 indicates the positive sense and physical 
meaning of all vector components of Table 2.1. Using eqns 
(2.18) in Table 2.1, eqn. (2,15) gets expanded to. 


m(U - YR + ¥Q) = mg^ + + 


T 

X X 


rQ(V + UR ~ ¥T) = mg + -Fa + 

^ y y 


m(W - UQ + VP) = mg^ + 


(2.19) 


To expand eqn. (2.17), rewrite its left hand side using 
the vector triple product, that is, 


/ £ X (ii. >< £ + tiL (il X £) ) ‘^^ 


V 


= / “(£•£) Pji dV - / r(r,w) dV + 


( 2 . 20 ) 


V 


V 


/ dV -/ r X r(^.^) dV 


V 


The first and second term in the right hand side would give, 
/ !i(£*ll) dV - / r(r.M) dV 

= (iP+jOfl-cR) / (x^+y^+z^) dV - 

V ^ 

# « • 

/ (ix+3y+kz)(xP+yQ+zR) P^^ dV 


V 



Table 2 
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Definitions of Sector Components 
of Equations (' 2 .'l5l and (2* fy')" 


iP. + 3P. + kP, 


(2.18a) 


for tiie aerodynamic force components. By a special orientation 
of the HZ axes these forces can be shown to be drag, side-force 
and lift respectively. 


iJ’m + jPm + kPn 


(2.18b) 


•for the thrust force components. 

for the components of gravitational acceleration, 


MOMENTS : 


il^ + 


(2.18c) 


(2.18d) 


for the aerodynamic moment components: aerodynamic rolling 
moment, aerodynamic pitching moment and aerodynamic yawing moment 
respectively . 


= iD„ + jM„ + kHn 


(2.18e) 


for the thrust moment components: thrust rolling moment, thrust 
pitching moment and tiirust yawing moment respectively. 


u =; iP + jQ + kR 


(2.18f) 


for the angular velocity components : roll rate, pitch rate and 
yaw rate respectively. 


= iU + + k¥ 

--P 


(2. 18g) 


for the linear velocity components; for-ward velocity, side 
velocity and downward velocity respectively. 


)ISTAITGE3 


£ s= ix + jy + kz 


(2.18h) 


for the components of the vector which locates mss particles 
dm inside XIZ, 



f igufe 2,2 i t iana of Vector Cumponcuts in the c ions o f 


Mot ion 
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* i{P / (y^+z^) P. d7 - Q/ 3jy p . d7 - R / xz 

^ JX ^ ^ 

3 {Q / (z^+z^) p dV - P/ yz P, d7 - R/ yz 

V* -tx ^ JX ^ 

k{R / (z^+y^) P. dV - P/ zz p, dY ~ Q/ zy 

-y " V V 

Phe integrals are now recognised as moments and products of 
inertia. Using tiieir common symbols, the right hand side of 
eqn. {2,21 ) becomes 

= - Qljy ~ 

These inertias can be calculated if the size and location of 
each mass in the aircraft is known. The third term in the 
right hand side of eqn. (2.20) similarly becomes 

/ £ x S!i(»-£) Pj^ 

V 

= / (iz+jy+kz) X (iP+jQ+kR)(Pz+Qy+Rz) P^ dY 

V 

= H PH+Iy^(E2-Q2).l^^ PQ+EQ(I,,-Iyy)) + 

Host of the airplanes posses a plane of symmetry along their 


p^ dY} + 


p^ dY) + (2.21) 
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longitudinal axis. If the X-axis is selected to lie in this 
plane of symmetry, then = 0. Using this property 

of aircraft and expressions (2.22) and (2.25), eqn, (2.17) 
in the component form would yield, 




Iyy« + = “a + “l 


(2.24a) 

(2.24h) 




(2.24o) 


Bqns. (2.19) and (2.24) form six differential equations with 
U,'V,¥,P,Q and H as six variables. 


2.5 ORIENTATION OP THE AIRCRAPT RELATIVE TO THE PHEL 
COORDINATE SYSTEM X* Y’ Z’ 

Consider system X' Y* Z' translated parallel to itself 
until its origin coincides with the centre of mass P of the 
airplane. This, translated system named X^ Y^ is shown in 
Pig . 2.5. 

The relative orientation of X Y Z to X^ Y^ Z-^ is described 
by means of three following consecutive rotations. 

i) Coordinate system 2^ is rotated about the Z^-axis 

over an angle i{) , positive as indicated in the figure. 

This yields coordinate system Xg Y^ Z^* The angle t is 
referred to as the heading (or yaw) angle. 
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ii) Coordinate system ^2 ^2 rotated about the T 2 ~axis 
over an angle 6 , positive as indicated in figure, 
yielding coordinate system X^ Z^. The angle 9 is 
referred to as the attit u de (or pitch) angle. 

iii) Coordinate system Xj Z^ is rotated about the X^~axis 
over an angle 9 , positive as indic 3 ,ted in figure. This 
yields coordinate system X Y Z. The angle 9 is 
referred to as the bank (or roll) angle. 

Yaw, pitch and roll angles are called Euler angles. 

2.4 COMPOIENTS 0? GRAVITATIOHAl FORCE 

Bqn. ( 2 . 19 ) indicates that there are three components of 

the gravit actional force mg. In the body -fixed system X Y Z, 

these components are mg^, mg„ and mg„ respectively, 

X y z 

It follows from Fig. 2.1 that, 

^'E = K^g = ig^ + jgy + kg^ (2.25) 

where the unit vectors K' and refer to the X* Y’ Z' and 
X^Y^ Z^ systems respectively. 

It is desirable to express g_, g and g_ in terms of g 

X y Zi 

and the Euler angles © and 9 (the yaw angle i}* does not 
affect components of gravity). 

From the relationship between the linear and angular 
velocities of 1 Y Z and X' Y* Z* systems, it follows that, 
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» -g sin e 

Sy ~ S sin 9 cos 9 (2.2o 

g^ ts g cos 9 cos © 

These relationships are now suhstituted in eon* (2,19) to 
obtain eq.n. (2.27) given below. 

2.5 HBVIBW OP THE EQUATIONS EBRI7BE 

The force and moment equations, eqns. (2.19) and (2.24) 
become, respectively. 


m(U - VR + WQ) -mg sin © + S' a + 

X X 

m(V + UR - WP) = mg sin f cos © + 
m(¥ - UQ + VP) = mg cos 9 cos © + P^ + 

Z 1& 

and 




'xx 


zz yy ■ 


lyyQ * 


I„R ~ + Uyy-Ij;x)f3 + 1X2® =11^ + % 


-ZZ xz 


(2,27) 


( 2 . 28 ) 


These eq\xations form six simultaneous nonlinear 
differential equations of first order in terms of 
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six variables** U,V,¥,P,Q and R. ' ^ ,/ 

Solutions of the above nonlinear equations are useful 
in certain applications, such as airplane simulators, 
accident investigation and the analysis of aircraft response 
to large distircbances. Most airplane motions of practical 
interest to the aircraft designer are steady state motions 
(that is, equilibriiun f3-ight) and perturbed state motions 
involving only small perturbations relative to some steady 
state motion. 

It is the perturbed state motion we are interested in 

r 

for our type of problem. Perturbed state flight is mathe- 
matically described by considering all motion variables 
as the sum of a steady state quantity (subscript l) and a 
perturbation qtaantity (lower case). Similar substitutions 
are carried out for the forces and moments, Phis is the 
process of linearizing the nonlinear equations describing 
aircraft motion about its steady state motion. 


Bqns. (2.27) and (2.28) also contain variables © and 9, which 
are related to P,Q and R by means of a set of three 
kinematic equations, which, in fact, should be invoked for 
complete trajectory of the airplane. But these are being 
excluded as they would not have a direct bearing on the 
problem of ^nthesis we intend taking up ahead. 
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Bqns. tz.jv) 

. ’ therefore, yield 

m{u - (y +v)(B ^^, 

t ^ (W,+W)(Q + )j 

- .ng eih(e + 6) + j ’ 

A f. + P 4- -P 

\ 

- if, 

“°"'^i+9)=oe(e,ie),p 


(2.29) 


^77^'^ ^ ^xx"^ Z2 ^ ^ ^1 +P ) ( +P ) 


2^ Z Z^ ^Z 

f '*I~ 

‘ ZZ-^yy ) (R^+r)(Q,j+q) = 

\ ^ h-^ ^T 

^^XZ ^ (^■,+P)^“(R,+r)^} = 

\ + + m (2.30) 


^ZZ^~'^XZ^'^^^yy'"^XX^ ^^1+P )(0 4. ^ T 

1 ^^(^1-^^)+Ixz(Qi+9)(R,+r) 


Define perturbations e' 


cos 6 


sin © -x. 


' xz''^r‘i>''‘-^l’^^>' = 

*"1, ■" “A + J’t, + “t 
9 to he such that. 


9 i. 1.0 


sin 9^9 


(2.31) 


This restricts the pitej^ 

roughly 15°, which is still • perturbations to 

View. "" sizable from a practical point of 


Expanding eqns. ( 2 . 29 ) ^nd (2 .n. 
ctions (2.31), it is fo applying restri- 

that the eqns . of motion can be 
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written as in Table 2.2. 
Remarks : 


i) Parts of tlie equations in this table underlined with one 
line are the steady state equations embedded in the 
perturbed state. Since, the steady state equations are 
inherently satisfied, they can be eliminated. 

ii) A few other terms are underlined with two lines. 

These are all nonlinear in nature, that is, they 

contain products of the pertixrbation variables u,v,w,p,q,r, 
© and 9 >. 

At this point, it is assumed that the perturbations are 
sufficiently small for products to be negligible with respect to 
the perturbations themselves. 


The majority of airplane dynamics problems are concerned 
with perturbations relative to a steady state for which : 


a) no initial side velocity exists : = 0 

b) no initial bank angle exists : 0^ = 0 (2.32) 

c) no initial angular velocit^ss exist : . 

!■, =«,=»!= = e, = = 0 

Introducing restrictions (2.32) and neglecting the 
nonlinear terms mentioned above, equations of perturbed motion 


simplify to ; 

m(u+ W^q) = ~mg cos 6^ + f^ 
m(v + U^r ~ ¥^p) = mg (p cos 
m(w - U^q) = -mg 9 sin + f^^ 


+ f 


T 


X 


+ f 


+ f„ 


’T 

z z 


+ f 



(2.35a) 

(2.33b) 

(2.33c) 



ti®as of Ifetiem with tfae ^strictions cose r cosji r 1 
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I p-I r=l. + 
xx^ xz A T 


lyyq = 


*" = 3a» + Urn 

zz xz-^ A T 


(2.54a) 

(2.54b) 

(2.54c) 


These eqiaations form the basis for most studies of 
airplane dynamics and control, 

¥e shall now proceed to determine expressions for the 
aerodynamic and thrust forces and moments to be substituted 
in the above equations- Though, experimental methods (using 
models) to evaluate these forces and moments have the 
advantage of achieving rather accurate values, in most 
preliminary design and parametric design studies theoretical/ 
empirical methods are employed. In so doing, an aircraft is 
considered to be broken into a number of components, namely, 
wing, body (fuselage), horizontal tail and vertical tail. 

The total force acting on the aircraft in some coordinate 
direction is then assumed to be the sum of forces acting on 
its components in the same coordinate direction. In estimating 
thrust forces and moments which depend on the power-plant 
characteristics, it is assumed that engine manufaetxurer* s 
data are available. 

2.6 STABIIITY AXES 

To derive and analyse aerodynamic and thrust forces and 
moments, it is convenient to use the so-called 'stability axes' 
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system of coordinates. For this, consider the airplane to 

be flying a symmetrical steady state straight line flight 

path. That means ^ = IL^ =7^ = 0 but and are not 

zero. The angle between the free stream velocity vector 

and is called , the steady state angle of attack. The 

stability axes X Y Z are obtained from the body axes X Y Z 

by rotating about Y = Y^ over an angle until X coincides 

with Yp (see Pig. 2.4). The new axes X_ Y Z are again 
^ s s s 

considered to be rigid ly attached to the airplane and thns 

move with it. It should be evident that the system X Y Z 

s s s 

is oriented relative to X Y Z differently for each steady 

state flight condition. Observe that in the X Y Z vsystem 

s s s 

U. = ^ Sind ¥. = Y. =; 0 which leads to further simpli~ 

' s -^1 ‘ s * s 

fication in the equations of motion. 

An aircraft in a steady state flight but with Y^ -1^ 0, 
is said to be sideslipping. (See Pig. 2.5). Sideslip angle § 
is defined in Pig. 2.5. The stability axes system in this 
case is defined in such a way that the X^ axis lies along the 
projection of the steady state velocity of the aircraft onto 
the 2L Z plane. 

2.7 PERTUHBED STATE PORCES AND MOMENTS 

Table 2.5 lists the perturbed state motion and control 
variables and the corresponding forces and moments. In this 
table 6^, 6^, and 5^^ are, respectively, the aileron, the 
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D,MB iti2a-£L£a^^ 



27 


elevator and the rudder deflections which are applied to move 
the respective control surfaces to cause desired motion of the 
aircraft. The question, which of the perturbed motion variables 
u,v,w,p,q,r and their time derivatives affect the perturbed 
aerodynamic forces and moments is answered by this table. 

This table in conjunction with Fig. 2.6 gives a systematic 
analysis of what takes place aerodynamically as perturba,tions 
relative to some steady state are introduced. It is clear 
from this table that it is possible to breakdown the perturbed 
state forces and moments into independent longitudinal and 
lateral directional sets. This, actually, is an offshoot of 
the fact that any inertia coupling is neglected, that is, any 
longitudinal motions which by definition lie in the plane of 
symmetry of the aircraft are independent of any lateral 
directional motions which lie outside this plane of symmetry. 
Henceforth, our approach/ study would be confined to only 
the lateral motions of the aircraft. We, therefore, branch off 
with eqns. (2.33b), (2.34a) and (2.34c) that comprise the 
lateral set and rewrite as : 


m(v + U^r) = mg cp cos + f^ + f. 


p- 

-I r 

= 1 a 

4* 

Im 


xz 

A 


T 

m 0 

r-I n 

= n* 

4“ 


zz 

XZ" 

A 


T 


(2.35) 


¥e are, therefore, left with v (or p fT*^ » P ^ 

1 

only variables of interest. Among time derivatives of these 



-Tahle Dependence of Perturbed Aei^dvnardc Forces and Moments on Motion, and Control VarlabUa 



1. All perturbations are considered relative to a sytaroetrlcal steady state: - Fj » 

^ a. »l«k» U ch« t.!.U indicate that thate 1. no effect to a fltat order of approalmaflon 




















*31/^ in stability 
^ axes 


\ 


“Z-Ht 


In vector notation; 

• U 4 -u + v + w, where V is the perturbed 
pi ^ P 

state total velocity, n is the steady state 


total velocity, while a,v,w are the pertuxhed 
velocities 


eC • arctan 


w 


U, 


1 } 

6 - nrctan ^ 


Interpretation of. Perturbed State, Velpcltifa .a^_Ag a | j £ i . 
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motion variables, experience has shown that v (or p) is the 
only one of any practical, importance. 

In addition are to be considered effects of perturbations 
in control surface angles. Here it is already assumed that 
longitudinal control surface perturbations (elevator etc.) 
do not affect the lateral directional forces and moments. So, 
we shall only confine to aileron and rudder perturbations 
6^ and 6^^ respectively. 

At this point it is also assumed that the perturbed forces' 
and moments depend only on the instantaneous values of the 
motion and control variables and not on the time-history of 
these variables (i.3., quasi-steady asstunption) . 


finally, it is assumed that it is reasonable to consider 
only first derivatives of forces and moments with respect to 
motion and control variables. 


With these assumptions, it is now possible to express the 
perturbed aerodynamic forces and moments as 


3F 


Ay 


3^ 


P + 


^^A 'V 

2L e + IX B + IX 

• P * ^ YV * 




3p 


ar ^ + 


3?, 

-.Jx 6 

36* ^A 


1"^ 

3P 
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31. 


31 


A 

5 


31 


3 1. . “ *>. 

+ ^r + ?57^A 

3p 


A 

A 
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A 


3H, 


“A = p ^ if ^ ^ 


SH, an. 3H, 


( 2 . 36 ) 
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It is observed that the j3, 6^^ and 6^ derivatives in 
(2.56) are dimensionless variables (i.e, radians). It is 
desirable to have all differentiations take place with respect 
to dimensionless variables. Division by an appropriate 
quantity converts these derivatives into their non<»dimensional 
form. Non-dimensional model of eqn. (2,36) would give 



where is the steady state forward velocity (along X) and 
b is the span (wing etc.). 
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Hewriting, in matrix form, eqns. (2.57) in terms of 
control and stability derivatives which are the quantities 
expressing the variation of forces and moments on aircraft 
caused by the disturbance of steady motion and are computable/ 
measurable, we obtain. 



forces Non-dimens iona,l 

and moments motion varia.bles 


(2.38) 

where q ( = t is the dynamic pressure and S is the 

reference (Wing etc.) surface area, 

* 

•j?or a complete discussion of control and stability 
derivatives, refer [15], pp 4,80 - 4.97f* 

Now for the lateral directional perturbed thrust 
forces and moments, it turns out that only the p variable has a 
significant effect on the airplane stability. We, therefore, 
introduce 

See Appendix A for definitions. 
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y 


(2.39) 


In terms of p-thrust derivatives, (Refer [15], PP 4.1)0), 
eqn. (2,39) is written as. 




sb 



(2.40) 


Adding eqns. (2.38) and (2,40) to form the total perturbed 
forces and moments of the aircraft, we have. 


f. +fn 


1^+1^ 


C C C G 
^pii ^r ^6, ^6, 


p r 6, 6^ 


0^ +C„ G C G G 

Uq Uq, n n Ue n^ 

qSb P ^ ‘^A ^R 



where, G . , G,., G . 2L 0 . 

yp J-p ^p 


(2.4T) 
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Combining eqn. (2.41) with eqns.of motion (2.35), the 
lateral directional small perturbation equations in the 
stability axes system (¥^ = 0) obtained are, 


m(v+U.r) = mgcp cos 9, + q.S(C f+C ■^+ 0^^ •“ + 


*6^) (2;42a) 


pb „ rb 

2b' °1 2U ‘^1^ ^ 

P P 1 r 1 0 * 6r 


^xxP~^xz^ = q^Sb(Cp -P+Cp l^+Gi ‘V^lc (2.42b) 

\ 


=«iS'=(0 -p+a -p+G 


.pb Q .rb_ 

~ />TT ~ 


P I, 


"p 




% 'h * °n. '^s’ 


(2.42c) 


2.8 STATE^VBC'TOR EQUATION Ob' MOTION 


* 

Substituting p (i.e., roll rate) = 9, r (i.e., yaw 
rate) = , v (i.e, perturbed side velocity) = U^ p and 

dividing eqn, (2.42a) by m, (2.42b) by I__ and (2.42c) by 
we get. 


« 

p ■ 

. is.. 

mU^ 

G 

•P • 

cr * 

~ 9 +• ’ 1 ' 

• 

= 0 


• • 

9 • 

qSb 
~ I 

^1 ' 

•p - 

' ^ ■'^A • 

-f^°l -^R 

X JL c JX 

= 0 

(2.43) 


XX 



XX 0^ 

XX 0^ 



• # 

9 . 

iSb 

I 

G - 
n,. 

'P ■ 

„ iSb . 

I ^n. A 

■ ?: 

= 0 



zz 

P 



zz 6 |^ 




where it is assumed that I^^, = 0 as both I^ and ^ ^ ^xz 
and that the total pitch angle 9=0 as the flight condition 
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does not involve any steep ascent/de scent of the aircraft. 
Terms involving the rotary derivatives** and the effect of 
control authorities on side force*** have "been omitted fojr 
simplicity . This does not pose any serious restriction Ofi 
the process of synthesis we intend discussing subsequently. 

The control and stability derivatives involved in eqn. 
(2.43) have been defined in Appendix A, 

Eqn. (2.43) finally represents a system of linear 
differential equations with constant coefficients that des- 
cribe the perturbed lateral motion of an aircraft. 


♦ • 

If ^2 *9, = 9 » and x^ =: , then the state 

of the system me.y be expressed as a four-component vector 
X, or 
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Similarly , the control deflections may be written as a 
two -component vector u, or 


U ss 


(2.45) 


The system (eqns. (2.43)) in the state-vector form may 


be written as 


X = X + B u 


(2.46) 


where 


5L§-. rj 

mu, 


I ^1 

XX p 
^Sb Q 

^zz 


(2.47) 


gSb „ 

I °1. 

XX 6, 


n 

XX o. 


(2.48) 


SlSfe Q 

I _ _ n c 
zz o , 


iSk ft 

I... 

zz Or 


The elements of the ^stem matrix A and the control distribution 
matrix B are constant parameters that depend on flight condition 
and airplane characteristics. 
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2.9 OUTPUT EQUATION 

Of the four elements of the state~vector defined above, 
we shall concentrate on <p and p responses of lateral control 
motion. It will be shown in Chapter III that application of 
either aileron or rudder by the pilot would result in a 
coupled response in terms of 9 and p, 7/e intend decoupling 
this response vector for his convenience. So, defining our 
output vector as. 



'^1 

72 


9 

z = 


L. 

our output description would 

X = c 

X 



where 

0 

0 

0 

C = 

1 

0 0 

0 


( 2 . 49 ) 


Its these matrices A,B and C obtained above that would 
be input to our algorithm of synthesis procedure and thus get 
useful results for the solution of decoupling problem. 


2.IO I'CEASUHBl'fflNT OP STATE YMIABEES 


In the sequel, decoupling is achieved by the application 
of a control law that involves state feedback. So, the state 
defined above must be sensed or' picked up' by some means to be 
fed to the error detector. 
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The roll angle 9 is measured by a vertical gyro which 
has its spin axis aligned with the airframe z-axis and is 
proportional to the voltage from the gyro pick off. 

Measurement of sideslip angle p involves measurement 
of the relative wind direction or the direction of relative 
motion of air as it passes over the airframe. This is usually 
accomplished by means of a vane, a probe, dual pressure pick- 
ups, or some similar device. 

The rolling and yawing velocities 9 and 9 are measured 
by the appropriate rate gyros by sensing the torque which is 
generated by the gyro due to an angular velocity input. 
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CHAPTER III 

/ by state feedback 

(.with reference to aircraft lateral dynamics) 

The purpose of this chapter is to state the decoupling 
problem and then, by proper state feedback, to find a .suitable 
control law that will not only remove coupling but will also 
enable the designer to assign system poles to achieve desired 
stability , Algorithm and computer program and the pioneering 
work done by Cilbert and Pivnichny to accomplish the above 
will be explained. It will be shown that aircraft lateral 
motions are coupled and need to be decoupled. Approach to 
decoupling problem formulation and solution will be general, 
i.e., applicable to multivariable nth order system. Reference 
will be made to two~input two-output fourth order aircraft 
lateral control system wherever desirable. 

5.1 heed to decouple 

All aircrafts have to be fitted with a control system 
that will enable the pilot to trim and maneouvre the aircraft 
in flight about each of its thre_e axes. On all conventional 
aircrafts, the controls are ’aerodynamic' ones; the moments 
required about the centre of gravity of the aircraft being 
produced as a result of changes in the aerodynamic forces 
acting on the various parts of the aircraft. In the majority 
of cases, the changes in the aerodynamic forces are produced by 
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means of flap-type control surfaces positioned at the 
extremities of the aircraft, the movement of the pilot’s control 
causing the flap surfaces to move relative to the main structure 
of the aircraft. There are usually three separate control 
systems and three sets of surfaces, one set for eontrolling the 
aircraft about each of the axes. ''The pitching control is pro- 
vided by the elevators, the rolling control by the ailerons 
and the yawing control by the rudder. The eleva.tors and 
rudder are usually part of the tail tuait and the ailerons are 
at the wing tips, in order that all these provide the longest 
possible lever arm about the O.G. of the aircraft, 

' It is generally desirable that each set of control surfaces 
should produce a moment only about the corresponding axis. 
However, in practice, moments are often produced, about the 
other axes as well., In other words, the effect produced by 
these is normally coupled. VJhereas the action of elevator is, 
by and large, independent of rudder and aileron, the la.tter 
two are interacting. For instance, the yawing moment produced 
by the aileron deflection and the banking moment caused by 
rudder action have generally been described as ’adverse- and 
may, at times, be viewed by the pilot as unpleasant handling 
qualities.^ The adverse aileron yaw also tends to reduce the 
rolling performance of the aircraft, because the yawing 
produces a sideslip and the dihedral effect produces a rolling 
moment opposing the desired aileron rolling moment. Similarly, 
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the rolling moment ineritabiy produced by rudder control can 
be termed -adverse' as it opposes the yaw effect that the 
rudder may have been put on to produce. 

The implication of the above discussion is that the 
coupling action of aileron and rudder controls not only would 
reduce their individual desired effects (and hence demanding 
the pilot to exert more) but would also require the pilot to 
use considerable skill in order to simultaneously manipulate 
these two inputs and successfully control the aircraft in 
lateral direction. 

Hence, to the extent that the control of the aircraft in 
the lateral direction gets easier for the pilot, enabling him 
to achieve better flight-path accuracy, there exists a need 
to decouple the cp and p responses. 

In, fact, many of the control and stability derivatives 
are an offshoot of coupling effects. Physical explanation of 
the aerodynamic mechanism of two of these relevant to our 
problem is given in Appendix B, 

Whereas the cross-coupling between the yaw control and 
the rolling moment and roll control and the yawing moment 
produced by the application of control surfaces can be avoided 
by a feedforward mechanism wherein parts of the total (i.e., 
gross) aileron and rudder deflections are cross-fed to rudder 
and aileron respectively, for the lateral directional decoupling 
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to b-o more effective, the aircraft motions in the roll and 
yaw modes generated as a result of the application of con- 
trols snould be decoupled. It is suggested that while the 
decoupling of moments can be achieved by neutralising the 
effect of off-diagonal elements of the open loop transfer 
function matrix, i.e., by so devising the cross-feed as to 
make the off-diagonal elements zero, the decoupling of cross 
rolling and yawing motions resulting in offsetting of yawing 
moment due to rolling motion (9) and rolling moment due to 
yawing motion (■/)) can be effected by diagonalising the 
system transfer function matrix by state feedback. The 
algorithm/ synthesis scheme made use of in the sequel attempts 
to achieve the latter. 

3.2 THE DBCOUPIIFG PROBLEM AID ITS SOlUTIOI [1,2,4] 

Consider that the time-invariant linear dynamical 
system S « {A,B,C} to be controlled has the form (dotted 
portion in Fig. 3.1) 

X » A X + B u 

3 ^ s C X (3.1) 

where u 6 ys ^ ® ^ real constant 

matrices of appropriate size. The transfer function of the 
system is 

H(s) » G(sl - A)""^ B 


(3.2) 
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a.ad tiie system output for zero initial condition is 
I(s) = H{s) U(s) 
or in expanded form as 


y^(s) -h^^(s) u^(s) + h^2(s) 0.2(3) + ... + h^j^{s) o.^(s) 

72(3) -h2.|(s) u^(s) + 1122(3) 02(3) + ... + Ji2jjj(s) %(s) 

( 3 . 5 ) 

* 

These equations are coupled , since each individual input 
influences all of the outputs. 

Def initioi^. 1 ; A system of the form ( 3 . 1 ) is decoupled (or 
noninteracting) if its transfer function matrix H(s) is 
diagonal and nonsingular, i.e,, 

y.j (s) = h^ ^ (s) (s) 

72(3) :s ^22(3) U 2 (s) 

( 3 . 4 ) 


Ve vshall nov examine some of the principal decoupling 
results used in this work. Necessary and sufficient conditions 
were developed by Falb and ¥olovitch [2,12] whereas further 
work leading to a complete synthesis of the decoupling problem 
was achieved by Grilbert [ 1 ] . 
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In what follows, the problem of decoupling will be 
discussed with, respect to the linear state feedback control 
law of the form 


U = f X + & V 


(3.5) 


where ? is an mx n real constant matrix, G- is an m x m real 
constant matrix and r is the new m input vector. ¥e will 
denote the above control law by the pair {P,G-) . G enables 
the inputs to be rearranged/reordered. The question of 
interest here is whether such a control law can be used to 
produce a closed-loop system which is both decoupled and 
stable. 


The desire to decouple raises several practical question, 
namely, (a) Is decoupling possible? (b) What is the class of 
control laws which decoupled? (c) What is the class of decoupled 
closed-loop systems? (d) What is the correspondence between 
elements of the classes mentioned in (b) and (c)? These four 
questions constitute the decoupling problem as it is treated 
in this thesis. 

Substituting (5.5) in (3.1)# we obtain 
X = (A + Bi') X + BG y 

The transfer function of the closed -loop system 
S(B,G} sCA+BB, BG, 0 ) , thus obtained, is given by 
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H(s,P,Cr) ss C(sl - A ~ BG (3.7) 

¥e shall give below the conditions on the open-loop 
transfer function, H(s), under which the closed— loop transfer 
function, H(s,F,G), can be made diagonal and nonsingular, so 
that decoupling is achieved - 

Define the non-negative integer d^ as 

d^ ^ min(the difference of the degree in s of the denominator 
and the numerator of each entry of the ith row of H(s))-1 

and the 1 x m constant row vector D^^ as 
d.+ l 

4 Lt S ^ H.(s) 

^ 0 00 ^ 

where Hj_(s) represents the ith row of H(s), 

Let 



Theorem 1 : The system (3.1) may be decoupled by using control 
law of the form (3.5) if and on]y if the m x m matrix D is non- 
singular. 

The above theorem signifies that whether a system can be 
decoupled or not is a property of its transfer function matrix. 
The dynamical-equation description comes into play a. role only 
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when the feedback gain matrix F is to be determined. Therefore, 
the controllability and observability of the system equation 
are immaterial here. 

Expanding (5.2) using leverrier' s algorithm and using the 
definition of d^ given above, we obtain 

A ^ B ^ O' (3.8) 

where denotes the ith row of C, 

Therefore, d^^ can also be defined from the open-loop 
system equations as the smallest integer such that 
G^A Bj^O'. In other words, integer d^^ = j where j is the 
largest integer from {1, ..., n-1 } such that B =0 

for k s= 0, , . . , j-i . 

If A^ B O' for all k 5. then d^^ = n-1. 

Similarly, we may define d^ and for the closed-loop 
system and obtain 

hi - C^(A + BP) ^ BG O' (3.9) 

It will be interesting to note that and are not affected 
by introducing the state feedback, i.e,, they are P-invariant. 

D efinition 2 : An P-invariant of S is any property of S(P,G) 
which for any fixed G does not depend on P. 

T heorem 2 : For any P and any nonsingular G, we have 
= d^ and D^^ » D^G. 
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y/i'fch. "fehe above preliminaries, Theorem 1 can be proved 
establishing its necessaiy and sufficient conditions as 
stated below : 

^.^QQQ.sit;;: — Suppose that there exists i* and Gr, such that, 
H(s,f',G) is diagonal and nonsingule.r. Then, J,, 


— 


- d. 

A ^ B 

d^ 



* 

n 

« 


A ^ B 

* 

« 

•d 

in 


a 

o 

— 


_ -J 


is a diagonal constant matrix and nonsingular. 

Sufficiency ; If the matrix D is nonsingular, then the 
system can be decoupled. 


If we now define another' matrix 


.«■ 


0, i 


Oj 1 


d ^+ 1 

V', 


m 


d +1 

m 


^then by confining our choice of P and G to 

P = a* and G = 


( 3 . 12 ) 


( 3 . 13 ) 


the system oan be decoupled and the transfer function of the 
decoupled system is given by 
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Or, equivalently 




C^(sI-A-BP)""'’ BG * diag(-j^) 

s ^ 

So, the above theorem specifies a control law which 
results in m independent subsystems, the ith being equivalent 
to dj^+1 levels of integra,tion, i.e,, 


yfCs) = Vj_(s) for i = 1,2, m 

s ^ 

v'“^' Such a system which is said tc be in Int‘=grator'-Decoupled 
(IJj) form, is not satisfactory because all the poles of the 
decoupled system are at the origin. So, additional feedback 
will have to be introduced so that while retaining decoupling, 
it is at the same time, possible to achieve the freedom to 
alter the pole locations (and hence the dynamic behavioirr) of 

* t-" 

the individual subsystems. To accomplish this, it is convenient 
to put the system into a canonical form. 
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3.2.1 Canonically Decoupled Systems 

Suppose, we perform a state transformation x = T x, where 
T is nonsingular, on the system S = {A,B,C} defined hy (3.1 ). 
The new system, thus obtained would he 

S = { I,B,G} = {TAT“‘', TB, GT~'*} . 

This motivates the following. 

Definition 3 : S and S are similar if there exists a non- 
singular n x n matrix T such that TA = AT, TB = B and C = CT. 

In other words, T carries S into 3. 

Remark 1 i If S and 3 are similar, their transfer function 

matrices are equal. Thus D = D and d^^ = d^^ for i = 1,2, m. 

Therefore, if S is .ID, S is ID. 

Remark 2 Simils-r systems are control law equivalent (CIS) 
implying thereby that a one-to-one correspondence can be 
established between and {?,&) such that 

ii(s,R,G) =H(s,S,G). 

The main theme of this subsection is to suggest that 
every ID system is similar to a canonically decoupled (CD) 
system. This means that if S = {A,B,C} can be decoupled, it 
is possible to find a CD system which is CliB to integral- 
decoupled S. Since, for a CD system (defined below) the decou- 
pling problem has a particularly simple form, it implies that 
the treatment of the decoupling problem for S is simplified. 
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D efiBiition 4 : S s{A,B,C} is canonically decoupled (CD) if 
the following conditions are satisfied : 


(i) The matrices A,B and 0 have the partitioned form : 


A =s 


A^ 0 

0 A, 


0 


G 


0 

0 


0 


A^ 


0 


B = 


0 


0 

0 

0 

0 


0 0 
0 0 


A 


u 


a: 


u 


m 


0 


A 


u 

■'m 


m 3214* I 


A 


u 


m+1 


0 


A. 


'm+2 


Ai Is 

4 Pm+1 - Pi, 


.U 


A. is p. X p. 
1 : 


m+2 


0 

b. 


0 

0 


•i-P 

0 


m 


m 


0 

0 


0 0 


m 


where p^ ^ i = 1, 


0 

0 

0 

• 

♦ 

* 

0 

m. ■ 


b . is p . X 1 , 


\ Pm+1* 1 » 


u 


•“m 


c^ is lx p^ 


^i ^ Pm+2 ^ 


(Pi 


0 ) 


Pm+ 1 ^m+ 1 


> 0, 


■m+2 


■m+2 


> 0 . I i • ■' 



Mob 
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(ii) Por i = 1, m the matrices A., h. and c. have the 

1' X X 


partitioned form : 


r 

0 I, 


0 0 


Ti is Ti * (d^+l). 


0 is r^ X r^ , 


» “ ! 


y. are diagonal elements 

f ^ 


of D = r ( say ) . 


0-^— fl 0 


where r. = p^-1-d. and is >. 0* 

^ iX. mLm, 


(iii) Por i = 1, m the systems == ^'^i'^i^^i^ 


are controllable, i.e., the p. column matrices 


Pi-1 

hi, Ahj^, A hj^ are linearly independent. 

m 


(iv) Let p = I Pi and t) =s •••» ^ 


partitioned n~row where is a p^^ row. Then, if 0 0 


and "0 is such that r) , . . . , r\ 

p-t- 1 P+P, 


are not all zero. 


then the 1 x m row matrix function ri(sI~A)'“'* B has at least 


two non-zero elements. 


Let <£ denote the n-dimensional space of n element 


row matrices. For i = 1, ..., 


m, define 
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X = {ri|Ti r]A^B^ =0 for k = 1 .,in, and j=0,. ., ,n-l} 

where is the ith column of B, 

lemma . j ; Assume S = {A,B,C} is IB and controllable. Then 
for i = 1 , m the following conditions are satisfied. 

(i) X, ^ is a row invariant subspace of A, i.e., 

T) G implies hA e X 

(ii) X. nX. = (0) for 3 = l,,,,m, j i / 

(iii) C.,C.A, ...» C.A are linearly independent elements ofoC.. 

XXX -1- 

The above lemma implies that there exists a linear ^ace 
'^m+l ^ such that the direct sum X ^ (+) -^2 ***‘^m+1 ~ 

let, jjj 

p. ^ dim X., i = m+1 and p = I p^ 

^ ^ i=1 

Therefore, p = n-p is uniquely defined by S although 
in+ 1 

X , . is not (unless = 0). Also, from part (iii) of the 

m-h 1 m-r 1 

lemma above it follows that p^ ^ d^^+l, i = 1,,.., m. 

¥e can now, state the following lemma. 

Lemma. 2 : Assume S = {A,B,C} is ID, contro .liable, and 

IQ+ 1 

p . 5 ^ 0. let T] e and write h e I 5-; » where 
m+1 ^ 

e X j, for i = 1 ,..., m+U IfC^^.^ ^ 0, then there exists 
at least two integers from the set {l,...,m} , say q and r, 
such tha.t T)A^B_ ^ 0 for at least one j S ^ and 

rjA^B^ i=- 0 for at least one 3 g ^0,,.., n-1 K 
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Replacing S by S in Definition 4 for the CD foim and 
using the results of two lemmas stated above, the following 
would stand proposed* 


Proposition 1 : Assime S is ID and controllable. Then S is 
similar to a CD system S, where ^ ~ s.nd 



In the proof of the above proposition, the matrix 




is used, where the rows of the p^ x n matrix fr^m a basis 
for Thus, the rows of Q (denoted by* say, q*, q^, •«.,!*) 

form basis for ok. . Q is used to carry S into S. To be more 
specific, we define 
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where q.* » • • • » q* are any linearly independent rows which 
i id, 

extend G^A, C^^A to form a "basis for 

Proposition t requires S to "be controlla"ble. To make it 
more general and so applicable to uncontrollable S7stem as well, 
we need the following lemma. 


Lemma 3 : Por the nth-order system S = {A,B,C} , let 

= dim ^ , where ^ is the subspace spanned by the mn 
columns j =0,..., n-1 ; k = 1 ,..., m. Then S is similar 

to S = {A,B,C> , where 


(i) 



A 


u 


A 


U 


B 



A IS n^ X n^ , 
c c 

a'^ is n^ X (n-n^) , 

A^ is (n-n ) * (n~n^) 
c c 

■nC * 

B IS n X m , 
c 


C = C^] 


C is m X n 


,u 




is m X (n~n ) , 

c 


(ii) = {a'^,B®,C°} is controllable, 

( iii) If S is II), S® is IB and r° = r where r is diagonal and 
nonsingular and represents the matrix L* 


The procedure required to form a CD representation of an 
ID system can now be summarised. Let S be an nth order ID 
system. Apply Lemma 3 obtaining subspace ^ and then a 
matrix T^ which transforms S into §. Apply Proposition 1 
replacing S by S®, which is ID and controllable. The matrix Q 
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in this case would be n x n . Define * 

^ c ]Pi 


T, 


Q 

0 


a+2 


0 


I- 

P 


n-n and 
c 


IlH-2 


Then direct ecalculatlon would chow that carriee 3 into the 
OD eyeten S, where p. » di» . 4, i = , Therefore, 

Tj I, would oari7 S into S resulting in the truth of the 
following theorem; 


Tia^ms-? .. . ! Ev-eiy M system is similar to a OD system. 
5 . 2.2 Principal Results 


In pursuance of the fact that the oa: uloal representa- 
tion of the system is more oonrenient to work with, the 
decoupling problem would be solved for the OD system and 
then by means of Theorem 3 , the results can be carried to the 
original system through back transformations. 

‘^-= - “ S is OD, the control law {p,s> decouples S 


if and only 

if 








0 

0 ... 

0 

0 




0 

®2 

. 0 ... 

0 

0 

'"2 


P = 

• 

• 

m 






• 

« 

4 

* 

» 




• 


* 

« 

« 

« 

f 


0 

0 

0 

0 

m 

0 

m 

[ 

where 

©i is 

1 X 

and 

JL 

is t 

* Pm+2" 

G 

= diag( 


X. 

i 

^ 0, 

i * 1 

1 ♦ 
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Suff iciencjr of this theorem follows hy direct substi- 
tution which shox-rs H( •,!',&) is diagonal. In fact, the ith 
dia,gonal element of H(s,f',l) is given by 

h.(s,f,(}) = c . (Ip .s-A^ -b.&. b-Y^.^,- (5.14) 

JL 11X11 11_L 

xhiere diagonal elements of matrin r stated in 

lemma 3 p8,rt (iii). 

'In consequence of Iheorem 4, we state the next following 
theorem ; 


'Theorem 5 ; Assiome 3 is Gi) and the control law {fjG} has the 
form indicated in Theorem 4. Then, 


(i) ±he ith diagonal element of H(s,7,Ct) is 

cc . ( s ) Y - ^ • 

h^(s,-,d) = » 4 = m (3.l4a) 


r^. r.-l 

where a.(s) = 3 "" - a -.3 

1 i1 


cc .. 




and 







a. =[ cr. 

X ^ xp. 



(ii) a^(s) = det(I^_ ,• 

(iii) 9. = ( o. -n;.. )!/'. , xirhere Y. is a p. ^ p. nonsingular 

wXm -J-* mXrn •JLi J- fju* 

matrj.x which depends only on and b^, and the 1 p^^ 

matrix tc^ = [*-■ ... 0 *** 

(ixr) The charact -ristic equation of the closed-loop system. 
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m 

q (s,P) where 

q:^(s) = det(I s -A,), i = m+1, m+2 , 

■^i 

Prom (3,26), h^(s, P,G-} may he interpreted as the 
transfer function of the system 3^(6^, X^) where 

= {A^, 1 )^, 0 ^ } is. a controllable, single-input, single- 
output system. So, single variable system theory, 
summarised in the following lemma, may be applied to 
prove the above theorem. 

lemma 4 : Assume S = {A,b,c} is single input, single 
output, order n and controllable, Then the transfer function 
of S(@, X ) has the form 


H(s,e,X) 




X 


where w(s) is a polynomial of degree n-t or less and 
’i'(s,o) = s -a.jS ^ - ... - o^.let OssCor^... ]and 

define the matrix K = [k^ ^n^’ where the columns 

4f 

kf = b, i = 1 n. Then K is nonsingular and 

©K = cr - -n;, i.e. , © s: (cr - 'n;)K'’^, 


Except for the form of (s), this lemma directly 


proves part (i) of the above theorem. Prom the form of A^^, 


it is clear that q-(s) 


d .+ 1 

ai(s) 


det(I s 

PT Pi-1 i 

«il® 


p . -r . 

“ 4 ^ s ^ where 
ir^ 


S QCjt t B j * S 4 o I 

we have used the notation of (ii). letting correspond 
to of lemma 4 verifies (iii). 


**These are given by leverrier ’ s expansion [l^J and cctopris© 
coefficients of the numerator polynomial of H(s,9,A), 
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We now estate the final theorem which embraces all the 
main results stated above,* 

Theorem 6 ; Assume S = {A,B,G } can be decoupled. If 
{F,G-} decouples S, the diagonal elements of H( *,?,&) have 
the form given in part (i) of Theorem 5 where the integers 
and r^ and the polynomials ct^Cs) are uniquely determined 
by S, and Yj^ = 1 f i = 1 » m. Furthermore, q(s,F) has 

the form given in part (iv) of Theorem 5, where a ,As) and 

iht 1 

“m+2^^^ are polynomials of degree p^_^^ and 

determined by S. The class of control laws which decouple 

S can be characterised by G- 6 and F G , where -§ is 

m 

an m~dimensional linear span and <? is a ( ^ + ^ip p) ■" 

i=l ^ 

dimensional linear manifold. More specifically, there 
exist matrices , . . , , i = 1, ..., m, and an 

(mp^_i_2) ~ dimensional linear space , which are uniquely 
determined by S, such that 

m 

G = I ^ G (3.15a) 

i=1 ^ ^ 

m p . 

P = -D-’a« + I (3.15b) 

i=1 k=1 

where F^ g ^ ~ •**' r^ , and = 0, 

Ic “ r 1, Pj[^* 

The form of H(*,F,G) follows from the CIS property 
between S and a CD system and Theorem 5. This would yield : 
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0 = d'U a * diag(X^, X^) (3.t6a) 



~~ {3.t6b) 

Results (3.15) are a direct consequences of (3.16), 
Theorem 4 and Theorem 5. Substitution of (5.1 6b) Into 
det(sI-A-BP) leads to the expression for q(s,P). 

The questions raised with respect to the decoupling 
problem get answered now. Part (i) of Theorem 5 tells the 
designer what closed-loop characteristics are attainable. 

Once h^(s,P,Gr) is specified by fixing cr^, eqns. (3.15) would 
give the corresponding control law. If the open-loop system 
is controllable, { P,(x } is uniquely determined (P^ would 
become zero). Part (iv) of Theorem 5 implies that closed- 
loop characteristic roots are all assignable iff = ^dh- 2"^ 

that would make =1. If "khe polynomials 

a^^^(s) and “124.2^^^ have any root whose real parts are non- 
negative, it is not possible to design a stable decoupled 
system. 

Theorem 6 establishes all the necessary data 
(p^, r^, ot^y D, Jj) needed for the design of 

decoupled multivariable system. Expressions, formulae and 
the general approach for obtaining this data have been given 
in the foregone discussion. 
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3*2*2 Computer Program : Description and Details 

The computer program accepts as input data the matrices 
A,B and 0 and produces as output all the data described 
above for synthesis. The main program, attached as 
Appendix G to this thesis, consists of 532 statements 
written in FORTRAN IV language and uses 17 standard SSP [20] 
subroutines. Names of these subroutines and their operation 
are annexed with the above appendix, A capability for 
holding problems of order n £ 25 and input-output dimension 
m £ 10 has been provided and so requires a large quantity of 
memory ( 100,000 bytes approx.). Since the original program 
was mm on IBM System 360, minor changes into it were 
incorporated in order to make it compatible with the ^stem 
DEC- 1090 . These modifications have also been quoted in 
the annexture. 

The progi^am after determining that decoupling is 
possible (det D 0) proceeds to transform the original 
system S into a new system 3 which is control law equivalent. 
Actually three such CIE transformations are made : 

. s S 4=== - S '^-==r' S 

Thus, a solution of the decoupling problem for S, a CD 
system which has a particularly simple and obvious form, 
signifies a solution for S. So, this solution when traced 
back through the same transformations wouM yield the desired 
results for synthesis. 
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The basic steps of the program are outlined in Fig. 3.2, 

SteB_.l and 2 ; Calculation of d^,D,A* 

Integers d^, ..*# d^ are determined by = j where j 
is the smallest integer such that ^ 0 for j = 0,..,,n--1 

and is the ith row of G. The B and A* follow from (3.tl) 
and {3.12), If det B = 0, decoupling is not possi.ble and the 
program ends. Else, it proceeds. 

Step 3 ; Conversion of S to S (IB Form) 

The transformation between S and S is given by 

: A = A - BB“'^ a*, B = BB”"^, G = C 

If f = BF + A*, C = BG, (3.17) 

a direct calculation would show that H(*,F,G) = H( •‘,F,G) 
proving, thereby, that S and S are CIE. A direct calculation 
also shows that d^ = d^ for i=l, ,..,m and 

B = I, A* = 0 (3.18) 

If a system satisfles(3. 18) , it is integrator decoupled since 

H{ s,F,G) = diag(-^;pj') 

B ^ 

step 4 : Controllable Form of IB System (S) 

The controllable part of IB system is separated from 
the me ont reliable part, if any, in this step. The contro- 
llability matrix H = [B,A B, 1%, B]is formed and it^ 

rank n^ determined. If n = n, S is alrc ly controllable and 
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Fig. tm% Flow diagram for main steps in decoupling pix)grani* 
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S is taken to be S. n, an n x n matrix I = [L^ I^] 

is formed where is a se t of linearly independent colijmns 
of H and L 2 "tke columns forming basis for the null space 

mmf 

Of H , i.e., they are orthogonal to the columns of . Then, 
with T^ = 1 “^ , the transformation 

A T^ = T^I, B = T^B, CT^ = C 

would bring the system S in the form given by lemma 3* Note 
that, 

P T^ = P, G = G (3.19) 

shows that S and S are CIP. 


Subroutine MPGR is used to find the rank, the independent 
columns and the orthogonal columns. 

Step 5 : Transformation to GB Form (S) 


This is done by first calculating Q. matrices from the 

d. 

portions of Form 0^, C?(A_) This forms 

part of Q . . The q 


O’. 

1 


qp are determined in the following 
i 


way, 


Form the controllability matrix H® for S°. From 
form the matrices i = 1, m where has the same 

columns as except for n„ col-umns i, i+m, . , , ,i+Cn„-1 )m 
which are zero columns. Then T)Hj^ = 0 is solved by sub- 
routine MPGR producing q . . . qp _ , which are chosen to be 

X ^ 

orthogonal to each other. Tlmt is, a basis for the( column) 
null space of the matrix 
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0 




1 


is obtained and the rows q. 


^i 


Pi-di-1 


for i =1, 

are obtained by talcing the transpose of these basis 
oolumns. 


m ^ f 


m 


Then 


,1 




is calculated, if necessary. 


Since, the rows of must extend the rows of 

to form a set of n linearly independent rows, the same are 

» 

fo-und by using subroutine MFGE with [Q^ ... 0^^^] . The rows 
of are the transpose of the resulting basis columns. 

Application of subroutine MPG-R also determines 

P^» P^+^ C= ^1^+^)- 

Tg is, therefore, structured and the transformation. 


: A Tg = T2A, I = ^ 

yields the GD form S, such that 


= F, ^ = a 


(3.20) 


Step 6 ; Determination of the S^n'^iiesis Formulae 

Since S is CIE to S, carrying {F,G> baeh through (3.20), 
(3*19) and (3.17) gives {F,G> for the decoupling of original 
system 3, 
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Form of and CD system (liieorem 5) stiggests that 

h. {s,P,G) correspond to applying control law ^6., to 


1' X 

Si = {Ai,bi,Oi> . For 8i = Vi(cyi-7ti), the itiUi,,... iti^^) 


and hence the coefficients in the characteristic poly- 

nomial for each are determined by using Devarrier* s 
algorithm. The E matrix (single variable theory) is formed 
and inverted to obtain 7. . When X. and © . are carried back 

i XX 

through the OIE between S and S, eqns. (5.15) and related 
data are obtained. The matrices in the print-out are 
given by expression of feedback matrix P of Iiemma 4. The 
matrices and are calculated as indicated on the flow- 
charts shown in Appendix C. 

5.3 AIRCRAFT MTERAI DYNAMICS: STATEMENT OF PROBIEM 


The linear dynamical equation for aircraft lateral 

motion was derived in Chapter II as, 

♦ 

X = A X + B u 
X = G X 

where u and x wei*© defined as two-component input and output 
vectors, respectively. 

The -transfer function matrix of this system is 


H(s) = C(sl - A)~^ B 
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where. 


sI-A = 




0 

qSb 

I 

XX p 


2Mk n 

S 


s 


-1 

0 

0 

s 


The determinant of the above matrix is 


A (s) =: 0 *8^4- (§ 

■^22 P 

= + PS^ + Qs 

= s(s+a2)(s+a^)(s+a^) 



(say) 

(say) 


It can be seen that 


C(sI-"A) ^ 


ZP« * ^ 1 


€L. SL§M Q .04. s Sb 
U, 1 s+ 

1 XX 22 Sj^ 


C„ -s" 


fTm. % 2 


•(S^- |§~ 0 -8+ ^ c )^ c + L 

1 p 22 p XX XX 22 


S- £Sfe Q .g. £Sk C . 

I ZX Ujj, iS2 Ug 


^(32 Q ^ 2 .§j 3 , rj 0 4,‘Xs^'Ul-- n -q 

1 P 22 P XX Oj, XX 22 p O 


"R 
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Sfi) 


Rs^+!Ps 


R's +T’s+K 

Therefore, 


?a^+¥s 
V s^+W* s+E’ 


(say) 


H(s) = 


Rs+T 

3 

s +-PS+Q 

Ys+W 

s^+Ps+Q 

1 

*^11 

^12 

R' s^+T‘s+K 

V* s^+W s+E* 



^22 

4 2 

s +Ps +Qs 

'"' 4 - 2 

sVps'^+Qs 



— 


where the values of P,Q,S,T,T,¥,H’ ,T* ,E,V' , aud E’ can he 
determined from any aircraft configuration pertaining to a 
fiying condition. 

So, inputs and outputs are coupled as indicated hy the 
transfer function matrix, i.e.. 


^ “ ^11 ^12 \ 

p = 1^21 ^22 

So, the statement of the decoupling problem pertaining 
to aircraft lateral dynamics is ; find a suitable control 
law {P,G-} which decoup les 9 and p. The same practical 
questions arise with respect to this problem ; 

i) is it possible to find a control law which decouples, 

ii) if decoupling is possible what freedom of choice for 

{f,G} exists, 

iii) what kind of closed-loop characteristics are attainable, 

iv) once an attainable closed-loop characteristics is 
specified, how is {F,Gr} determined. 
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All these questions will be answered once the data 
derivable from matrices A,B and C pertaining to a particular 
aircraft and flying condition is determined by running the 
above explained computer program. 


This would be taken up in the next chapter. 



CHAPTER 17 


APPMGATIOH OP BECOUPLIHC THEORY ^ 

(Achievement of Stable Decoupled D7namics) 

This chapter aims at illustrating the applies cion 
of the decoupling theory explained previously in arriving 
at numerical solution of the lateral dynamics decoupling 
problem as applied to a high-performance aircraft f=imilar 
to X-15 that was produced by North American Companj? as a 
research airplane for NASA. The aircraft configura'^ion 
and aerodynamic characteristics under stated flying 
conditions as taken from a NASA Report [6] and relevant to 
oiir problem are given in Table 4.1. Sy incoiporating the 
values of parameters from the computer output into the 
relations stated in the last chapter, complete synthesis of 
the decoupling problem, raised therein, is obtained. Step 
inputs of aileron and rudder deflections are simulated and 
the respective tp and p time responses in both open— loop 
and closed-loop cases have been drawn. 

4.1 COMPUTER INPUTS 

The A and B matrices, computed from (2.47), (2.48) and 
the data of Table 4.1, are 

-0.0176076 0.005367 

• 0 0 

^ = 2,6719776 0 

4.1260132 0 


0.3490658** -1 

1 0 

0 0 

0 0 


See foot-note next page. 



71 


fable 4»i Aerodynamic Characteristics of an M?ariced 
Airplane Similar to Z~15 


b,ft 

g,ft/sec^ 

I^,slug 

I^a^slug 

m, slugs 
q,lb/sq ft 
S,sq ft 
U^,ft/ sec 
Altitude, ft 
Mach number 

0-, ,per radian 

P 

C ,per radian 

Up 

C ,per radian 

• 

0, , per radian 

^A 

(3, - , per radian 


k • • # « 


n • * « ^ 


> 0 -0 m 


! ,per radian 

n6 


« 


>^6 ' 
R 


per radian 


22.36 

32.2 

5,021 

67,159 

390.4 

200 

200 

6,000 

125,000 

6 

20 

0,015 

0.31 

- 1.0 

-0.075 

-0.15 

0.08 

-0.108 
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B = 


0 

0 

•13.359888 

1.0647776 


0 

0 

-26.719776 
- 1.4374497 


The matrix C was already fixed in (2,49) as 



1 

0 


0 

0 



integers N (=4) and M (=2) are fed in as the sizes of the 
state and the input (and output) vectors respectively. 


4.2 SrUTHBSIS 

The computer print -out (attached as Appendix D) obtained 
after running the program for the design of decoupled system 
explained previously is referred. In connection with deter- 
mination of a control law {P,Cr} which decouples, the four 
questions raised in Sec. 3.3 are now answered. 


i) Since matrix I), uniquely determined by the dynamics of 
the lateral motion, i.e., A,B,C ,. is found to be non- 
singular, necessary condition of decoupling is satisfied. 
Hence, it is possible to find a control law that would 
decouple cp and p responses. 

This element of matrix A while, deriving the equations of 
motion had yielded zero by neglecting the nonlinear term wp 
in eqn, (ii) of Table 2.2. Its inclusion gives rise to the 
teim a in A whose value under the stated flying conditions 
is appreciable and hence may not be neglected. Though a, 
primarily concerns longitudinal motion, some inertial 
coupling is normally found in case of high-performance air- 
crafts.. 
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ii) The computer output gives 


h-SUB-t, 

i*e. , 

d, =1 

D-SIIB-2 , 


11 

P-SUB-1 , 

i*e. , 

CM 

II 

P-SIIB-2, 

X « @ 

P 2 = 2 


The fact that P-SlIB-3 does not appear in the computer 
output signifies ths.t the completely unohservahle subsystem 
found in the structure of canonically decoupled form is not 
present. This, in fact, would happen in most of the cases. 
Thus, p^ = 0. 

Therefore, r^ indices would get. 


Ii 

0 

using = Pj_-1-dj_ 

rj =0 


rj = 0 

(By def., * r^ = p^) 

and r^ = 0 

(By def., * * the open- loop sys. is 
* 


G ontro liable ) . 

Because, = 0 for i 

= 1, ..., (m+2), i.e., 4, by Theorem 5 


part (i) of Chapter III, the polynomials, a^^Cs) = 1 for 

m — 1, ..., 4". 


The class of {P,Gr} which decouples as given by relations 
(3.15) would be 


0.16556 

- 0.12020 


'0.56069 ^ 


G = ^ 1^1 ^ 2^2 


0.28035 ^ 
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and 


(P^ is not defined as system is controllable). 


j^. 23266 - 0.00683 

■0.00005 - 0.00397 <y^2 

+ 5.59530 ^22 

+ 0.00208 ^^21 

1.21633 + 0.00342 a,. 

{ ■ ^ 1 . 

0.00003 + 0.00288 0^2 

1 ~ 2.79770 022 

- O.OO 1 O 4 02^ 

-0.00044 - 0.12046 ^ 

0.00987 - 0.38800 

+0.13544 ^21 

0.00022 + 0.08746 

- 0.00494 + 0.19400 

- 0.06772 021 

^2’ ^11’ ^12* ^21’ ^22 

being arbitraiy real 


numbers denote the freedom of choice with P and G. This 
freedom enables the designer to place the closed-loop poles, 
theoretically, anywhere in the complex s~plane. 


iii) Elements of the diagonal matrix H(s,P,G) from (5- 14a) 
of Theorem 5 are. 


hj,(s,P,G) _ ^ 


i = 1 ,2 


where the polynomials and ^2 would be 




s -ar^^s-cy^2» 


^^2 ~ ^ ”^21 ^”^22 
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Therefore, 


H(s,P,G) = 




0 


0 


s 3-022 


(4.1) 


With i*s and cr’s to be chosen arbitrarily, the above 
•tells the designer what closed-loop characteristics are 
attainable. Once these are selected, a particular set of 
h^(a,P,G) is specified and the corresponding control law 
{F,G} determined. This follows. 


The open-loop poles** of the unaugmented aircraft in 
question whose lateral dynamics are coupled and unstable are 
found to be 


= 0,00449 

■•r 

and = -0.01105 + j 1.78698 

The desired stability is represented by the roots of the 
controlled system which were chosen to satisfy the military 
handling qualities specifications [18], These were selected 

V 

as follows ; 

= -0.0546 
= -0.695 

and = -0.346 + o 5.14 

Given by the eigenvalues of A, these are calculated by 
using IMSL subroutine EIGRP. 


t 
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Comparison with, roots of the basic (tmaugmented) 
airplane shows : (a) that the lateral motions of the airplane 
are now completely stable; and (b) that the damping ratio 
of the complex pole pair is much higher for the controlled 
system. 


Since, the roll motion of the aircraft is non- 
oscillatory by nature whereas the yaw motion is oscillatory, 
the simple poles selected above will be attributed to 9 
response while the complex poles will be associated with the 
p response. The selection of roots and their appropriate 
assignment to the 9 or p response also bears testimony to 
the fact that while, on one hand, the resulting ratio 

= r 3 + 6 ' :o ' 34 di: 3+0:693) Ts^aVcB+t " ) > 

with a % O'- and b>>a, which is recommended by design 

considerations, a time period of 2 secs, given by the complex 
roots above, on the other hand, is considered fairly accep- 
table for the transients of p-response. 


With desired stability, our system would, therefore, 
be described by 


H(s,J?'',G) = 


N 

( s+0. 0346 y(s+b. 693 ) 
0 


0 


( s+ 0.546- j3 . 14) { s+0.346-y 14) 


(4.2) 
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which, hy comparing with (4.1), yields 


=- 0.7276 

cr ^2 = “ 0.0239778 

( 4 . 3 ) 

= - 0.692 

0-22 = - 9.979316 


11 

II 

M 

V' 

arbitrarily 

( 4 . 4 ) 


then, in terms of value parameters in (4.3) and (4.4), the 
control law fP,Gr} which results in a stable decoupled lateral 
response gets fixed up and takes the form, 


P 


3.3548126 - 0.0019437 

- 1.5779528 0.0018741 


and 


0.0145368 - 0.5781974 

0.0199584 0.1890622 


G 


0 . 1 6556 
- 0.12020 


- 0.56069 

0.28055 


A similar set of gains can be computed for any set of desired 
roots. In a complete analysis, the roots may have to be 
chosen selectively to insure reasonable gains. 


The system matrix of the closed-loop system is found to b 


A+BP 


’- 0.01761 0.00537 

0..0 0.0 

5 0.014601 - 0.0241079 

I 9.9663755 - 0.0047635 


0.34907 -1.0 I 
1.0 0,0 : 
-0.7273205 0.0009735 I 

- 0.0132108 - 0.6748597 I 


Its eigenvalues and hence the roots of the characteristic 
equation are found to be 


0 . 0346226 , - 0 . 6921571 , 


- 0.34620 + 33.138997 
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which tally with the chosen roots- This corroborates oxh* 
synthesis of decoupling and stabilising controllers- 

4,3 SIMUMIIOIT 


Having obtained mathematical models of the system, in 
both Open- loop and closed-loop cases, we are now in a position 
to obtain their time responses- by simulating any specific 
input signal (say a unit-step). 


Referring to Sec, 3,3 and Table 4,1, the open-loop 
ti^sfer function matrix of the original system, after 
computation of various coefficients involved therein, would 
be given by. 


1.064783^ - 0,0717s 



-13.35988s^ + 0.23524s 
+57.96814 


-1.43745s^ - 0,1434s 

26,71978s^ + 110.24614s 
- 3.84081 


where u( s) = s( s-0, 00449 ) ( s+0,0l 105-3 1 .78698) ( s+0.0l 1 05+3 1 .78696) 


If a unit step input of aieleron deflection (i.e, , 

= i) X now siffiula.ted on the system, the resulting cp and .p 

M, 3 

responses would be given by 


„ _ 1.06478s - 0,0717 

^ ■“ s( s-0.00449)Cs+0. 011 05-3 1.78i^98)( 8+0.01105+3^.78^98) 


and 

g ^ -13.35988s^ + 0.23524s + 57.96814 

{ 3-0. 00449 ) ( s-fO. 01 1 05- j 1 , 78698 ) ( s+0. 01 1 0S-fj 1 . 78698 ) 
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If, on the other hand, the system is excited by a unit- 
step of rudder reflection (i.e, = ^), the corresponding 

Jx S 

9 and p responses would become. 


-1.437453 - 0.1434 

s( s-0 . 00449 ) ( s+0. 01 1 O 5 - 3 1 . 78698 ) { s+0. 0l 1 65 + 3 1 .78698' 


p ^ 26.7l978s^ + 110.24614s - 3.84081 

s^{ s-0. 00449 ) ( 3+0. 0i 1 05 - 3 1 .78698) ( s+0. 0l 1 05 + 3 ] .78698) 


iProm (4.2), the transfer function matrix of the finally 
obtained closed- loop system is given by 



For the unit-step aileron deflection, the output response 
of the system is 

cp = 

And, the output response for the unit-step rudder 
deflection would be, 

p = 

A complete time response expression for each of the above 
responses in frequency-domain can then be obtained by splitting 
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them into partial fractions and then carrying out laplace 
transform inversion mder zero initial conditions. This 
was achieved by writing a computer program. 

# 

Time-responses drawn have been shown in the adjoining 
figures. The open-loop responses (Fig, 4.1 to Pig. 4.4) are 
oscillatory by nature and indicate that the stability of 
lateral dynamics in the case of unaugmented aircraft is 
unsatisfactory. The decoupled 9 and ^ responses are shown 
in Pigs. 4.5 and 4.6 respectively and have been made 
completely stable. Whereas the closed-loop p response is 
found to have a settling time (t„) of 4.5 secs, approx., 
the 9 response turns out to be slightly sluggish, which, 
if desired, can be improved by manipulating the dominant 
root. 



31 
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FS9.4.4. Step response fsr^CinpulS^^J 




9 Cinpu.l<5A) 
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OHAPTER 7 

COMGLUSIORS 

For the aircraft lateral control problem considered, 
Gilbert’s comprehensive theory for the decoupling of multi- 
variable systems by state feedback provides a constructive 
procedure to determine the classes of feedback and inter- 
connect control matrices that give a decoupled structure. 

His algorithm, translated into machine language, is readily 
adaptable to any p rob lem. It is felt that this scheme would 
provide a better synthesis procedure than other methods 
because of the relative ease of solution of the pole place- 
ment problem. 

In arriving at the solution of the problem, our treatment 
does not take into account the transfer function due to 
electro-hydraulic power devices (comprising of servo-valve, 
hydraulic jack, motor etc) demanding deflection of control 
surfaces. As, in practice, this forms part of the control 
loop, its inclusion, that changes dynamics of the system 
and hence its system matrix, would pose no problem when 
obtaining the solution because of the generality of the 
method. 

While deriving the final form of equations of motion, 

control derivatives 0^^ and 0^ were neglected for simplicity 

P ^ 

which, in fact, is slightly unusual. Much as we would have 
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liked to retain these terms, their numerical values did not 
exist in the aircraft data.. But their inclusion would only- 
replace one of the zero elements in the system matrix hy a 
finite value and, therefore, would, in no way, prevent us in 
achieving synthesis- 

Whereas in oxir case, the free elements (cr's) of the 
feedback matrix (F) were fixed by specifying closed loop 
poles pertaining to desirable military handling qualities, 
the problem can be tackled alternately C4bj in a way that 
would optimise these free parameters with respect to a spe- 
cific criterion like minimising a cost functional involving 
error and control amplitudes. This amounts to chosing 
by solving an optimisation problem for S^» 

The effect of parameter variations and disturbance inputs 
on the aircraft system have not been considered as the method 
used makes no provision to take them into account. In pra- 
ctice, the system, which is decoupled for nominal parameter 
values, may not remain decoupled when parameters undergo 
variations. However, if the variation are not large, the 
system may be assmed to behave in an' approximately decoupled' 


fashion. 
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APPENDIX A. 

CONTROL AND STABILITY DERIVATIVES 


The control and stability derivatives involved in 
eqn. (2.58) are defined as follows. Their usual dimensions 
/^re deg**^ or rad”^ . 


Symbol 




3 P 

3C-, 


30 . 


R 


C. 


P 


= 36 ^ 
30 

:l 

“ 3 P 


0 

30 

^ i S L 

’"a 

^ng 3p 

P 


n. 


30 . 

■38 


R 


R 


30 

n Q 


Definition 


Variation of rolling moment coefficient with 
sideslip angle (i.e. dihedral ^angl| ) . 

Variation of rollii^ moment coefficient with 
aileron angle (i.e. lateral control power). 

Variation of rolling moment coefficient with 
rudder angle. 

Variation of side force coefficient with 
sideslip angle. 

Variation of side force coefficient with 
rudder angle. 


Variation of side force coefficient with 
aileron angle. 


Variation of yawing moment coefficient with 
sideslip angle. 


Variation of yawing moment coefficient with 
rudder angle. 

Variation of yawing moment coefficient with 
aileron angle. 
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Symbol 


Definition 



3 0. 








30 


n 




30 




30 . 






30 


n 





Variation of side-force coefficient with 
roll rate. 


Variation of rolling moment coefficient with 
roll rate. 


Variation of yawing moment coefficient with 
roll rate. 


Variation of side— force coefficient with 
yaw rate. 


Variation of rolling moment coefficient with 
yaw rate. 


Variation of yawing moment coefficient with 
yaw rate. 


Variation of thrust induced yawing moment 
coefficient (due to normal force) with 
sideslip angle. 




9 (|^) 


Variation of side-force coefficient with 
rate of change of sideslip angle. 


Variation of rolling moment coefficient with 
rate of change of sideslip angle. 


Variation of yawing moment coefficient with 
rate of •change of sideslip angle. 
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APPENDIX B 


CONTROL DERIVATIVES 

n 


AND ,£ 




•'R 




a) Yawing Moment due to Aileron (C ) 

“■6 


Consider a positive aileron deflection (Pig. B.1)j on 
the i^tarhoarcb wing, where the aileron is down, the lift of 
the wing is increased, and so are the profile and induced 
drag. On the aileron -up wing, however, the lift and hence 
the induced drag component is reduced. Whereas + and - AI 
would cause a hanking moment (as desired) that would star- 
hoard the aircraft, the difference in the drag acting on the 
two wings produces a yawing moment that would turn the air- 
craft to port side. Por this reason it is described as an 
adverse aileron yaw. 


If, on the other hand, the ailerons are being used to 
turn the aircraft to port, then the resulting yawing moment 
acts to starboard it as shown in Pig. B.2, 


Values of 0 are difficult to estimate and it is 

generally best to use experimental data. Its magnitude can 
vary widely from one airplane configuration to another. 


b) 


Rollii^ Moment due to Rudder 



) : 


A physical explanation of the aerodynamic mechanism by 
which rolling moment due to directional control is obtained is 
presented in Pig. B.3. 
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When the rudder is deflected, the resulting aerodynamic 
force 'P. ’ on the vertical fin and rudder produces a re- 

My 

quired yawing moment about the C.G-, of the aircraft through 

the distance arm X™. along the X - axis. Usually, this 

Vg s 

resultant force acts at a distance (along the Z^-axis) 
above the longitudinal axis through the C.G. and consequently 
produces a rolling moment as well. Because, the purpose of 
the rudder is directional control, this rolling moment due to 
rudder must be seen as an undesirable side effect. 


To develop an expression for C, , associate the force 

due to rudder on the vertical tail with a ’ side-force 

derivative* , C__ , so that 

R 

y( rudder) ^ 

Therefore, the rolling moment about the X (stability) axis is 

3 


^(rudder) " 

y s 


s 


R 


R 


( non-dim ensionalised ) 


so that; 


"R 


= C 


Uy 


0 is normally known is terms of vertical tail area and 
lift curve slope. 




APIEITDIX 


n 

'j 


C? MAI'T A” 


1 5 

i [) 

XI 

X? 

X 5 


/ .'*> 

-1 


5 ^ 


,1 ^ 

^ i'j 

j ^ 

4 r* 

4 1 
4 X. 
■'j3 

"I i 



itr siD'^ 
! a '*»* S 1 i'J ; 
>.; HFO 
rl D I li 

■•t 

i)U‘ 


> i( ! I 4 f W ■ I ,* f * » I I I J J Ji ’ J » ^ * -f f /■ f * ' S • Hh iU ■8,«J' ^ 

el.uilu’ Xoi7ifv( 10,25) 

.«; 'S f ^.>1 ;.b ^^' ( 25 , 25 3 , H ( 25 , 1 0 J , CB -\xt 1 1 !i , 25 1 
i<. : i> ) u ■■•. A J’ r. 1. ( 2 5 , 2 5 3 , B n I J C 2 5 f 1 U ) , c: f .1 L ( I < i , 2 5 ) , I. H f Hm ( 2 7 5 ) < 
Kxl ^ ^ 25 3 , Khb ( 25 , 25 5 , T1 ( 26 , 25 3 , IPC i 2 3 

, iSVtdH '=^l2,3X,2M‘'!=,'t 2,3Xr3.5HZLKO TDr.ERArjC£s,FB. 1 3 

iFf4,Fo.n)GU TO 999 
Prd.MT 6, 4, 2, tot, 

-JRhATC 6 1 FI 1,5,2X3 3 

J 20 1=1, N 

‘I % j '. A j'teJC'fr V \ 1 4 «. (• 


t •.!« 


F ... 

nu 


nu 20 1=1, N 
HEAD*, CACI,J3 
f’O 21 1 = 1, .4 
Rfi.AU*,C8(I,J3 ,J=UM) 

130 22 1 = 1, F 
READ*, (C(Ifd3, 0=1, s'U 
F 3R.nATC10X,BHA MATRIX! 
■'^T.,T 26 

. i:5,'CAC.I,J3,0=l,!->) 
kATHIX) 

1 77 

■J 1 = 1,'- 

1 h, (HCI ,03 , J = 1 ,•*<) 
t i--'- >'■: UOXfiOC OATKIX) 

■ ■■.i .T 29 
„ . . > ■“ -1 /, 

, (C(! ,03 ,J = l,f!3 


iVTi'iT 
R.) 
rvT 


25 

I 


ii 

'1 i 


AI'RAV 
AkRAY 
A r\ K A Y 


( 2 » ‘ 

(2, 

(2,; 


Ii i‘ 


f ' R T ' I 
f = 1 
CAl.l, 

CAJ.u 
C.Ai.'L 
0=0 
2 ai;.,. 

CAT.L 
CAf.Li 
,< = ■,;; 

!< = h', + 1 

IF (RtK)-T':)U3 
JF («()<) + roil 3 '(3, 
(•'.-lO 48,60,40 


■ iCFY 
iiCPY 
G.-.FT’O 


A, 25, 25, A, A3 
'•! , 2 6 , 1 0 , b , B 1 
.0,10;26,C,C 

-03 


( b , . 

(C,.l ,CI,M.0,03 


V. , i , V. .1. , ' " , « 

(Cl ,AO,R,l 


,'■*,03 


I? . . 

I OC I) 
!'b 4 3 
5 ' ( 1 , R 3 
IF ('l- 
T=T + l 


iO 
J 
5 = 1 . 

= RCK3 

13 48,60,45 


41,41,43 
43,42,42 


>5) 

, ICTLil 25 3 



!• 

^ 


H 7 


hU 


7 1 
d 7 


6 9 

957 

9 


y, } 

ll I 
^^7 


?1 


/ 

/ j 

/ -4- 

f 'i 
-y.} 

v'i \i ) 

^ 7 i 


) (A| A0| ATenp^^r 
CnTK^iP.ACi, 

'.i 7 t.'- Fku (Cl f AJjrPfl f 

C !-l-JJ 4P,43j9 
I r s (^' 


'I , M ) 


••■■n T'i 35 

i'Ki.'i 1, *4*^ 

f * '■.•.■, A t t14ij PRUGKAM KRROR) 

(■ ,1 f n 1 0 
C-U.l,, G’-PRD (A, 

C'^Ru ''C?V 

J-J + 1 

C ‘■' f 1 i ,1 !.j '■ F K u (Cl 

IF C 

I' 5F 

K.irF'li' (2M' ■) MATRIX IS SINGULAR) 

,I()‘'”FS‘V=1 
r,;i 'J.\i f)S 

CAf u AFkFi' J0,D,n) 

CAiai “CRv (i:-,i!rRV',N-,f<,o5 

C'\hu "I^vv 

JUT£SI=U 

(DET-TUb) 61,bl,65 
(DEr + T.lb) 65,62,62 
V,.. TO 57 
PHidT b7,L>FT 

FJk,AAT C2SH nETERMIftjAbT DP D "ATRIXsEI 0 ,3 ) 
CALL ARKAV C 1 , M , 4 , 1 0 , J 0 , D , D ) 

CALL " “ i*' V f 1 - M vi 1 n < ,”» A T ... u ri j fti \( ■> 


TP 
( F 
Gil 


L ARKAV Cl ,M, 4, 10, JO, D,n) 

L hRKAY (l,'%^,l0,10,DlAiV,in'^V) 
'■'.Al' (10X,6Hr.J '.Ai’PIX) 
ir 68 


FuR 

PRlir 68 


968 I = i , 'i 
KC-'if 15, (0(T,J),..J = l,h) 

lAT (l,0X,loHD IVVERSP MATRIX) 

HI a dR 

R 6 9 i s 1 , i‘i 

. !’ 1 5, (OX ’ ai,j) . jsj ,M) 

■./■■r (76 0-5bh-,l2, JH=,I2) 

;,n2 ' 


fM) 


F 3K 


!' T 


•nl 

1 


-50 h 

t 

'^4 1 a’ Rii. , Carol K) 

48,70,550 

rs ! 

7*.i.,L i^C»Y fC,I,CI ,''aFa>) 

CAi.'u 4r‘Pkii 

ra-b 'tC?y (K'i,ci,i,u,o) 


T = J + 1. 
i.F 


, CJ-IUC,L)-U 7J,72, 
7 3 F=1, 


48 


V 


St A^Hl ,KJ=CHr) 

1'- ) 7 !i , y 8 0 , 4 H 


rrc i‘ 

T = I + 1 

nJ 'iM 7U 

lJ(a,J3HA SI’AR 'MATRIX) 

PhL;) 80 
I ' i i y H 1 T s 1 j, rl 

('AT iT t5,(ASTARll,d),j3l,M3 

C.Mh ARRAV (2,iaRr1.0,25,A5TAR,ASTAR) 

CALL ARRAY (2 , 6 ,iV , 10 , J 0 ,IUNV ,DINV ) 

(H,DI.'"V,BBAR,M,M,M) 

(n6AH,A6IAH.Te?lP,N,' 

r K n a k n aj ai i 


Mi 


f u 


CILL 

CALL 


GFFKl) 

G'''‘pP*l> 




.NV 



C'iLb 'iCPy U) 

CAta., , AsrAft/i'E'ip,? 

cat.l -- '■" -- - 

C\hU 
CAT,!., 

CAf,,U 


*' I 

;■< ?. 
ii a 




. O J , M « * I, L - , , , . , 

t J ,^■,.V,25,?5,TEMP,TPMP) 
tU v,.'.,,25,2t), ABAft,rtaAR) 
,25/ 3 0,aBAk,BdAa) 
,'-',P,10,25,C8A(?.CaAS) 


h 

o 5 


><0 

)? 


A P R A y 
AB RA it 

ARRAY (3 , . 

APhAY (1 10,25, C8AR,. 

PijK'iAr (1 ox, .1 2 hA EUR f ATRIX) 

Pal. Vi' Kj 
1-3 982 1 = 1,0 

E81UT 15. CrtHABt r,J3 ,J=3 ,^i J 
KJK^.A'J’ CfiiX,J2HB BAR fUTRIX) 
rPTOT H2 
oa 93]li l = t,H 

E>k 1 .VE 1 5 , ( R8 AR ( T , J 1 , J=1 , 
format C3uX,3 2HC BAP '''ATHIX) 

Phi or K,? 

aJ 84 Isl.p 

PRiaT 15, ICUARCI, J) ,J=1 ,43 

FORMA r (lC’X,22HD-TA'VeraiE4A5TAR MATRIX) 

PRI.jf 85 

DO 86 I=1,M 

PRINT 1 5 , ( fERP ( 1 , J 3 , J = I , 'O 


CABCULArE T1 MATRIX 


t I I 

i .a 

I T 

J (SO 

1 2 ! 


K=n 

CALIi 

CAf.B AnRAV 

Cam, 

CAMj 
C41,t! 

1 . j I'M 
i=avi 

cao t'Ai'r 

t ; V - ( 


ARRAY ( 2 ,> ,F, 25 ,in,«HAK, 8 BAR 3 
hah A V (2 , << , 0 , 25 , 25 , A 8 AR, AI 5 AR 3 
' CRY CH-Ul ,AD,f.,{!,U 3 

aClit. CUfi'i ,0Y!‘.-!,F,0 ) 

a 'H'f'. t Ai’,'TRFP,a,.M 

1 M=!,- 


('1 h>aMa,H,K+j,.M, 
“ ' '08,109,48 

a , 0 3 


.M,a,0,H»M3 


.... ... . 1,3*/M^ 

» = t +'■ 

Ci'aa, mcpy (All, ATE:Mp,N,a,o) 

C Ai.u (r-piO) tAbAR,ATKMP,an,N,4,M3 

•■a la 100 

C aa, ''FOR (-MT ♦M,M,TQB,4C,IR0W,tC 
IF taC-M 122,320,48 
I 0 ( 2 3 =a 

CAhi, '-CPy C/U4AK,ATIL,A,f;,0) 
cai.a, -'iCPY (bi:.AR,HTlIi,r,M,0 3 

r ii 1 i . a c? « V f . y. . i'j . 1 n . O'; . r ft B 


r % f n J- J. *J g i* f ’■* w J 

cai.a, -'iCPY (bi:.AR,HTlIi,r,M,0 3 

(2,y',0, 3.0,25,CaAR,CBAR 

:iUR,CTIT..,F,4,a) 

T1 ,0,0,4,03 


C.AJ.O ARRAY (i 
C«».L ".CPV tC 
CXjub SC0.?i 

r.Aia, ucb‘ 


] 21 

1 2 ? 
J 2.1 


SC0.?i (Tl,0,a,4,o) 

UCb\ (Tl,1.0,N,0) 

.•■CPY (T1 ,Fl,b,M,O,03 
T (3S!! KYKTEa-i IS CnMTRObbABJ 
1 *? 1 


Cal or r;j 

'■a 3'U 190 

^ 12 3,:JC 

■ (310 SYSTEM IS Ni 
■ (El l,,0,.X,M,0) 


* nj ,1 

Pkll''.r 

F'U-.rtM'i' . 

Call SCO a 


rWT 


i 4 ATRIX=I) 

4 C=,I 23 



•n 


62 


9!] 3 
t<3 


65 


B6 


CALL MCPy {C,CPAR,M,ri,U) 

CALL LMPRD(DIhV,ASTAR,TEMP,H,M,N) 
CALL APRAy C1,M,A.,25 j5,TEMP,!r|f'Mp) 
CALL rtPRAi' (1 r^,^,25,2S,A6AR,ALAR) 
U»A.,l,25;10,86AR;buARJ 
C A L L A P R A y ( 1 , M , L , 1 0 , 2 5 , C 8 A P , C B A R ) 
NJh.^Ar C1UX,12HA bar f-^ATRIX) 

PRI.'^r HJ 
nj 982 1=1, R 

PRINT 18, (A8ARtI,J3 ,J=1 ,iM3 
RJHMAT (l0X,12HB BAR HaIrIX) 

PRINT 82 
DU 98 3 1 = 1, fi 

PRINT 15, (BBAH(I,J),J=,l,M) 

FORMAT CiUX,12HC BAR MATRIX) 

PRINT B3 
DO 84 1=1, M 


PRINT i5,(CBAR(I,J),J=:l,N) 
FORMAT^ (1 OX, 22HD-INVEPSE + AS 


PRINT 85 
DO 86 1=1, M 

PRINT 15,(rEKPtI,J3,J=l,N) 


STAR MATRIX) 


c calculate ti matrix 

= 


I '»■! 

1 .>l 
!. :t8 

1 09 
uo 


121 


1 22 
1 23 


CALL ARRAy ( 2 , 0 , M, 25 , 1 0 ,BBAP ,RBAR} 

CALL ARRAY C 2 , L , N, 25 , 25, ASAR, ABAR) 

CALL MCpy (BBAP,A6,N,fi,0) 

CALL SCLn (R,0,N4M,N,L) 

CAr.L C'^iTPA ( An,Tr>(P,U,il) 

I.J 101 
.1 = 00 

CML KAUn (TEMP, J,H,K+J4M,N,0,N>t<M) 

IK lK-(/',-l)*.‘i) 108,109,48 
". = 1 + '' 

call LCPY (At), ATEMP,N,M-0) 

CU,L iJi-'PRD t AbAR,ATEMP,AO,N,N,M) 

U.l TO 100 

CALL KFGR ,T0L,NC»IR0W,IC0L) 

IF (NC-L) 122,120,48 

I prM+2)=o 

CALL MCPy (AHAH,ATIL,N,M,03 
CALL MCPy CLLAR,BTIL,f>,M,0) 

CALL ARRAY ( 2 , H , N , 10 , 25 ,CBAR,CBAR) 

CALL F>CPy (CBAR,CTIL,M,N,0) 

CALL SCLA ('n,0,N,N,6) 

CAT,L UCLA (T1,1,0,N,0) 
call >.cpy (ti,fll,n,n,o) 

FUHMAT (350 SYSTEM IS CnUTROLLABLE Tl MATRIX=I) 
PHI or 121 
r,0 TO 199 
PRINT 123, NC 

format (3iH SYSTEM IS NUT CONTROLLABLE NC=,I2) 
CALL SCLA (ELL,0,N,N,0) 

CALL AHkAY Cl,M»M,fl,275,25,H,H) 

PRINT 125 



1 ‘i? 
Id'.' 
1 D I 
ib?. 
,15S 


FORMAT ClOXfHHA TILDA MATRIX) 


199 


2on 

9 2 ') 1 
y 2 ■' ) 2 

i, 

1 


V. 1 .) 

1 1 

<i 1 >. 

2 13 

2 1 *1 
21 / 


DO 159 i = l,fj 
PRIRT 15,CATrL(I,J),J=l 
PRT.^iT 160 

FORMAT {10X,14HB TILDA 
(>{) 161 1 = 1, « 

PR1:^IT 15. (fifIL(l,J),J=l 
PHI^T 162 

FORF.AT C10K,14HC TILDA 
DiJ 163 1 = 1, M 
PHIOT 15, CCTIL(I,J),J = i 
CALL AHRAT (2,«,N,25,25 
CALL ARRAV C2,L,H,25,10 
CALL ARkAX (2, M, 0,10,25 
CALL ARRAY (2,0,0,25,25 

CALClilATF 0 MATRIX 


MATRIX) 

,M) 

MATRIX) 

,ATIL,Arili) 
,ariL,BTlL) 
,CTIL,CfIL) 
, n, ri) 


HCLA (H,U, NC»‘M + irCI) + l,NC,0) 


219 
2 2 !) 


1 = 0 
iPfe;=o 

ISi)M=0 

CALL XCPy (ATIL,AC,1 ,1,WC,WC,>^,N,0) 

CALL XCPY CBTIL,BC,1,1,.VC,M,R,M,0) 

CALL XCPY CCTIL,CC,1,1,M,RC,M,0,0) 

CALL SCLA CO,0,N,N,0) 

00 200 K=1,A 

IS!iK=IDCK) + l+ISUM 
IF (ISUM-AT) 201,250,9201 
PRILT 9202 

FURi^iAl' (IHH The SUM ( D-SU8*I + 1 ) FOR 1 = 1, M IS GRFATER THAN NO 

CO 1') 4R 

r = r-H 

= '' 

CAi,,L ..Cpy ()iC,AU,NC,M,f>) 

CAL.t; LMRA (AO,TewP,NC,M) 

21 t JU = 1,«'' 

1= LJ 

O J.; L ^ i < A D t ( I e on > V J , ri , K + J , M , N C , 0 , N c * M + 1 D c I ) + 1) 

if'' (K- i'.C^ri) 212,213,48 
CALL ICPY (AO,ATeMp,N-M,02 
C\LL GMPRI) CAC,ATEMP,Ao,NC,NC,M) 
b'l ID 210 
« •J = Ci'.’C-l)4M 
D.i 214 L=0,Nri,M 

CALL SRMA Ch,0,NC»M+ID(I)+1 ,NC,I+L,0) 

CALL XCPY kTIL,CC,l,l,M,NC,M,N,0) 

CALL RCPY CCC,I,CI,M,LC,03 
[l=ID(I)+l 

CALL^RAfe^tcf ,1 ,H,NC*H+J,l,NC;OfNC*M+IO(I)+l) 

CALL RADD CClrl,a»IPE+J/i»NC,6,NC) 

CALL MCPY CCI,FN,1,NC,05^ 

CALL GriPRD (EN,AC,CI,i,NC,NC)„^, ^ 

CALL mFGR (H,NC*M+IDCI)+l ,NC,T0L,JR,IR0«,IC0L) 
IPCI)=3C-JR+1DCI)+1 



, ^ PKifrr 9220,1, iP(i) 

92^ > FUP,«'r (7H P-SUa-,I2,lH=,3X,r2) 

IPE=lPEflP(l) 

IF CIP(I)"ID(I)-1) 48,230,221 
/'2t CALL ARRAY ( 1 , NC^t'M+lOC t ) + l , NC , 275 , 25 , H, H) 

PHIMT 9221,1 

9221 FORiviAT (25H SOLUTIOiSi OF ETA^Hi FOR l=,I23 
X9C=OC*-44ia(l3tl 

DJ 9222 L=1,IRC 

9222 PRIRT 1 5 , CHCL, LI ) , Ll=l , NC ) 

FRI.^T 9127,CtCDLCL),L=l,tsl) 

Mi^NsNC-JR 

DO 223 J=1,^JNN 

CALL SCLA Cfi^,0,l,NC,0) 

EUICOLtJR+J)) = l 
DO 222 L=l,JR 

222 E.niCOLCL) )=H(L,aR+J) 

223 CALL HAUD Ca\ , 1 , Q , ID( I ) +1+J+IPE-IPC I ) , I , «C , C 

230 IF CI-M) 231,270,48 

231 ISUM=U 

DO 232 JslIIjM 

232 ISUMslSUH+IDCJ)+l 

IF (IPE+ISUrt-NC) 201,250,48 

250 laltl 
XPCI)aIO(I)+l 
PRIRT 251,r.IPCI) 

251 FaH.v!Af (7ri P-SI!B-,I2, 1H = , 12) 

CALL RCPY (CC,.1,CI,M,NC,0) 

1 £=jn(l}+l 
DO 252 f,=i,I£ 

CALL KAOn CCI,1,Q,XPE+L,1,MC,0,NC) 

CAbO <CHY CCI,F«, l.NCfO) 

/5,> CAt.L 0-1PHU CK?J,AC,CI,i,MC,NC} 

rv’(-:=jeF+T»-'C 1) 

XF (l-.i) 250,270,48 
IF aPF-KCJ 271.275,48 
CALL OCPV £g,T6^1P,^C,^C,0) 

ChLu fiFCR CT£PF,«C,NC,r0L,JRAMK,IRa’4,IC0L) 
XFCJRA jK-IPE) 48,272,48 
CALL ARRAY ( 1 , RC ,A!C , 25 , 25 , TEMP , TEMP) 

P R I v X 9 27 2 

FUR MAT (i6H SOLD r I ON OF Q»X) 

DO 9273 Lal,NC 

P8IRI 15, (XE^P(L,l.l),Ll = l,NC) 

PRI;^»T 9127,(ICOL(L) ,L=1 ,<'^) 

A OryiMsHC-XPE 
DO 274 

CALL SCLA (E('4,0, 1,RC,0) 

f'ruicoLCXPe+j5)=i 

130 27 3 L=1,1PE 

273 E.'J(ICnL(L))=TEHP(L,lPEfJ) ^ ^ 

274 CALL RADD t £0 , 1 , 0, IPE+J , 1 ,NC ,0 ,NC ) 
lPCM + t3 = !MC-IPE 

-oMMsM + I 

PRINT 9274,MMM,IP(M+U 
9274 FORMAT 17H P-SU8-, 12, 1H»,I2) 


27-1 

271 


9 27 2 
9 27 3 


,NC) 


2V> IF C'iC-fJ) 27b,2eO,48 

CALL SCLA {TFMt’,0,L,N,0) 

flCCl“>'iC + i ’ r r r 

CALL ARHftY C 1 , L , N, 25 , 25 , TE«P,TFMP> 
00 9276 Js^'^CClrN 
92/6 l£LPCJ,J)=l,0 

CAIJjAHRaY (2j-w,N,25,25,TEMP,TFMP) 

uQ 277 JJ=t,:-!C 

CALL SCLA (fe;N,0,l,‘^rO) 

J~d J 

CALL RCPY (Q, JfEH,NC,NC#0) 

27/ CALL RADD CEM, 1 ,TEMP, J, J , M ,0 ,N ) 

GO XO 281 

28) CALL .ACPV (g,TtMP,N,M,0) 

281 CALL MCPY ( tEMP , T2 , 0 ) 

CALL (r£i^P,N, DETER, Cl, EN) 

CALL ARRAY ( I , N , N, 25, 25 , T2 , T2 ) 
PRINT28? 

282 FORMAT C10X,9HT2 MATRIX) 

DO 283 1*1, N 

283 PRINT 15,(T2(I,J),J=1,N) 
IFCDETEK-TOL) 2§4,284,287 

284 IF(OETER+TOL) 287,285,285 

285 PRINT 286 

206 FORMAT C22H T2 MATRIX tS SINGULAR) 

GO TO 48 

287 CALL ARRAY C2 , L .N , 25 , 25 ,T2 , T2 ) , 

CALL GMPRO CCTIL, TEMP, chat, M,N.N) 
CALL Gs^PRO {T2,BT1L,8HAT,N,N,!V!) 
CALL G^PRO (T2,ATIL,AHAT,N,N,fO 
CALL Kei (AhAT,ELL,N,N,0) 

CALL G“:PPO ( ELL , TEMP , AHAT ,N , N.N ) 
CALU ARRAIE U , L , L , 25 , 25 , TEMP, TEMP ) 
call AnHAY ( 1 , A , N , 25 , 25 , AHAT , AHAT ) 
call array Cl ,f4,M,25, 10,BHAT,8HAr) 
CALL AKRAY (1 ,M,N,10,25,CHAT,CHAT) 
P K I 1' 2 8 R 

IHA P-jiH-AT (10X,nHT2 INVERSE MATRIX) 

C,0 289 T = 1,N 

^^9 PRI -iT 15, tTEMP(I,J) ,J=1 ,N) 

PRINT 290 ^ ^ 

29',) FORMAT C10X,12PA-HAT MATRIX) 

00 291 I=1,N 

291 PRINT 15,(AHATCI,J),J=1,M) 
print 292 

292 FURMAT C 1 OX, r2HB-HAT MATRIX) 

DO 293 1*1, N • 

291 PRINT 15,(BhATtI,J),J=l,M) 

PRINT 294. 

291 FQPmAI (10X,12HC-HAT MATRIX) 

LO 295 1=1,0 

295 PRInT 15, (CHAT(I,J) ,v) = 1 ,N) 

m 

Z CALCOLATE the G sub 1 MATRICES 

400 DO 406 II=i,M 

CALL SCLA (G,0,M,M,0) 


b i}U 

5Dl 


f = 1 1 

CftLL CAOD 
PHINT 405,1 

•liih format C6H U-SU&,12,2X,btlMATRlX) 

CALL ARRAlf ( 1 , , M , 10 , iO ,G,G) 

1)0 406 J=1#M 

406 PRIf4T 15, (GCvJ,K3 ,K=1,M) 

I CAiACnLATE THE PI X K CONSTANTS AND J MATRICES 

CALL ARRAX (2 , ^ , N , 25 , 25 , AHAT , AHAT) 

CALL AKHAY ( 2 , M , M , 25 , 1 0 ,6HAT , HH AT ) 

rpF=o 

CALL XCPX CAHA’l,AI,IPE+l,IPE + t,IPCl),IPCn,N 
CALL XCPY (BHAl,BI,lPEtl,I.IPa),l,N,rt.U) 
CALL SCLA (BI,0,IPa3 ,pU) ,0) 

CALL OCLA CRl,l,0,IPtI),OJ 
CALL DCPY (Al,EN,ipCI3,0) 

SUMsO 

IPI=IPCl) 

DO 502 K=1,IPI 

SUMssSUH+ERCK) 

QIKdIsSTJM 

J»0 

CALL SCLA 

CALL GRPRU CPi # BI , SN . IP ( i ) , XPC I ) , 1 ) 

CALL CAOD ,KK,I^CI)-J,lP(I),l,0,IP(n) 

call'" SCLA flENiPl ,0,N,^.03 _ 

CAlb .'CLA (TEMPI, OIK(J).IPCI 1,0) ^ _ 

r A L G ' ? K n ( H I e A 1 , 0 , 1 P ( 1 5 , 1 P ( 1 3 , 1 P (1 3 ) 


502 

50 J 


5 ‘ ' > 


505 


5 1 


pU 

bl 3 
05 U 

9514 
514 
5 I 
529 
521 

523 

524 


Ai,f» G 'pH n c H I { A I , y , it; I l R 1 1 J , if' ( 

41, Li ijr'Slit-jf U , 'lEMPl ,Rt,IPll3 ,iP(133 
Ip- 1 504,510,46 

r xr.L G ■i?Ri)fRT.,Al,0,IP(n,IP(l3,IP(I3 3 
r \l,ij uCPV (.(.», EM, IF (13 ,0) 

t'oj.i. — ;» 

i,)J 5o5 K=l,lPi 
SUM=SUM + EfjCK3 
vIK(J + 13 = SU vco + n 

c ALL ARR A Y C1,IPCI3,TP(I3,25,25,RI,P.I3 
:)13 9513 K=l,IPI 
ILI 9513 L=i,iei 
SUM=RICR,L3 . 

IF (SUM-TUL3 513,513,514 
IF (SUM+TOL3 514,9513,9513 
CD''TlMUt: 

(d TO 520 
I i I y 1 S - T 

KaRMt C'2Dii TRACE METHOD FAILS FOR I=,12 3 
PRIM’I’ 52l,I,IP(l) , 

FUHMhT (Bii I = ,I2,3X,8nP-SlIB-I=,I23 

PHIHT 524,1,K,0IK(K3 , _ 

FORMAT (2tl y , 12 , 12 , 1 Hs , El 0 , 3 3 
GO ru 600 


0,03 





9 57 4 
525 

5 2n 
577 

b 37 

9530 

9532 

9531 

53J 


bi^ 

5 33 
5 3'> 

33o 

537 
5 4 1 . 

550 

56(1 


ciiK'Xi-aJt: 

CALL tCt^Y CFK,VaPCn / rp(l),0) 

Caf,L ARRAY (1,TPCI),IP(I),2S,25,FK,FK) 
RRIrT 052a 

F’mMAT tl<)X,15HTHE K MATRIX IS) 

OU 525 K=1,1FI 

PRIRT 15. CFK(K,L)-L=1,IPI) 

CALL ^^IiyV CV,lP(I),Dt:T.EM,EM) 

IF (DhT-TOLJ 526,526,530 
IF (DFT+TUL) 530,527,527 
PRliir 520 

Fa((.‘lAT (21 H K MATRIX IS SUGULAR) 

(TO I'll 4ti 


KaO 

FOPRAt (tuX,15HTHe V MATRIX IS) 

PRIOT 9530 

CALL AKRAI Cl.IPCI) ,IP(I) ,25,25, V,V) 
DO 9532 !,= 1,IPI 


PRIRT 15.(V(L#Ln,Ll=:l,rPI) 

CALL ARRAY C 2 ^ IP( I ) , IP ( I ) r 25 , 25 , V, V) 


K""K + 1 


CALL SCLA 
CALL HCPy 
CALL SCLA 
00 531 Lai 
eM(IPE+L)3| 
CALL RADD 
CALL GHPRD 
CALL GMPRO 
CALL G0p.7i) 
CALL AH K 'll 
CALL AHHAY 
pHr-.l 63S, 
PHT'.T 53) 


(EMI ,0,M,N,0) 
(V,K.E/»,fPCl),TP(I),0) 
(EM|0, 1,25,0) 

£’R{L) 

(eM,l,EMl,I,l .9,0,M) 
(l>I^'V,EMi,EJl ,M,M,!^) 
(FJI,T2.TEMP1,M,R,,.0 
(reRPi,Ti,Eji,M.N,N) 
(1 ,i**,i'i,iO,25,EJl ,EJ1) 
(1 ,m,.n, 10,25,EMI ,EMl) 


(150 MATRIX H a) 

K, K is a i, P ( X ) + 1 •• H 

r Kl M 1 5 35 , KKhA 
I'M 53 3 Lal...i 

P R I T 1 5 , ( i; R 1 1 L , L 1 ) , L I a 1 , 9 ) 
FOP OAT Cl OX, 12) 

PRIM!' 535,1 

format (150 KATPIX J a) 

PHTw'l 5 36 
XkKK=IP( I )+l-K 
RfilwT 535,FK;<K 

00 537 Lai,« 

PRIPT 15,CEJ1 (l.,LU,Llal #9) 
IF (K-IPCD) 9531,540,48 
IF (1-3) 541.550,48 

1 ?Fa1 PF + TPd) 


T=I + 1 

GU TO 50 1 

lPE;alPt>lP(I) 

IF (IPE-NC) 560,10,48 
1 = 1 + 1 
GO TO 503 





c I 1 
o2C 
n/l 

i'2'i 

^2i 

b 0 1 

603 

604 
60*5 

b06 

9606 

607 

608 

610 

999 


IPJ = UHn 

IF CJ-0) 677,627,626 

THlsIPl 

c;u TO 620 

IF (IHI) 18,623,621 
lFCA66(yiKfIHI5 )-TUL3 622,622,625 
TRr = ItU-l 
GU TO 620 

IF 9525,9525,10 

PC(I«I+1)=1.0 

OJ 601 IA=1,IR1 

PC{IA)=-OIK(IRI+l-IA3 

CALL PRaOCPC,IRI+l,PO,PE,POL,IR,IKR) 

FORMAT t22H SUBROUTINE PROD FAILS) 

PRINT 6U5,I 

pur,v!atc41h NON-zEpn eigenvalues of matrix AHAT-SUB-,12) 

IF (lER) 606,606,610 
PRIM 9 606 

FDRf'iAT (lOH RE;al PART, 3X, I 4HIM AGIMART PART) 

DO 607 IA=I,IR 

PRINT 608,PO(IA) ,PECIA) 

format (ei0.3,3X,E10.3) 

IF(I-M)9525,9525, 10 
PRINT 603 

IF CI-M) 9525,9525,10 
END 






' 
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Caleulate matrices 



ilote ! The iiumlserg on the flowchart correspond to statement 
numbers in the program listing. 













107 


SUBROUTINES PACKAGE 


The names and functions of the subroutines called by the 
main program listed above are appended below. 


1. ARRAY : 

2. M3PY : 

3. ROPY : 

4. GMPRh ; 

5. Miinr : 

6. SChA : 

7. GMTRA : 

8. RADD : 

9. MPGR : 

10. PC LA : 

1 1 . XCPT ; 

12. SRMA ; 

13. GMSUB : 

14. CADD ; 

15. DOPY : 

16 . PROD : 


Convert data ARRAY from single to double 
dimension or vice versa. This subroutine 
IS used to li^ the user program which has 

^^ension arrays and the SSP subroutines 
which operate on arrays of data in a vector 
fashion. 

Copy entire matrix. 

Copy specified row of a matrix into a vector. 
Multiply two general matrices. 

Invert a matrix. 

Set each element of a matrix equal to a given 
scalar. 


Transpose a general matrix. 

Add row of one matrix to row of another matrix, 

Por a given matrix(M by R) E (a) determine rank 
and Itoearly independent rows and columns (basis) 

(b) find an orthogonal basis for the null space 

of B. 

Set each diagonal element of a matrix equal to a 
scalar. 


Copy a portion of a matrix. 

Multiply row of a matrix by h scalar and add to 
another row of the same matrix. 

Subtract one general matrix from another to form 
resultant matrix. 

Add column of one matrix to column of another 
matrix. 

Copy diagonal elements of a matrix into a vector. 

Calculate all real and complex roots of a given 
polynomial with real coefficient. 
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MODI FIG ATI OUS IN THE PROGRAM 


1 , All the DIMENSION statements were brought on the top 
of the listing as per the requirement of DE0-1O9O* 

2. READ statements were made foimat-free for the sake of 
convenience and to avoid confusion, 

3, Device numbers (Data Sets) from the READ and WRITS 
statements were removed. 

4. Since, the program itself did not have any halt state- 
ment, an additional statement 

IF (N.EQ.O) GO TO 999 

was inserted. Execution is terminated when the Input 
file reads 0 for N, 

5, Hollerith field specification was modified in certain 
formiat statements for printing of headings and variable 
names, 

6. Do loop index modifications •(to differentiate from the 
subroutine parameters) were performed to remove 
warnings , 



aEFBRBlTGBS 


1 1 1 


Papers 


1. E.G-. GiToerti 'The Decoupling of Multirariable ^stems 

by State Tariable Feedback', SIAfi J. COITHOI, Vol. 7, 

Feb. 1969, pp. 50-63. 

2- P.l. Falb and ¥.A. ¥olovitch: 'Decoupling in the Design of 
Multivariable Control Systems' , IEEE Trans, on Auto. 

Control, Vol. i!C-l2, Dec. 1967, pp. 651-659. 

3. E.G. Gilbert and J.R. Pivnichny: 'A Computer Program for 
the Synthesis of Decoupled Multivariable Feedback Systems' , 
IEEE Trans, on Auto. Control, Vol. AC-14, Dec. 1969, 

pp. 652-659. 

4. J.R. Pivniohny: (a) A Computer Program for the Design of 
Multivariable Decoupled Feedback Systems (b) Design of 
Decoupled Multivariable Control Systems, Systems Products 
Division, IBM Corporation, Bndicott, Eew York. 

5. Eugene Cliff and lutae; 'Application of Geometric Decoupling 
Theory' to Synthesis of Aircraft Da^teral Control System' , 
Journal of Aircraft, Vol, 9, Bov. l 972. 

6. W.E. Oehman and J.H. Suddath: 'State-Vector Control Applied 
to Lateral Stability of High Performance Aircraft' , FASi 

TN D-2894, 1966. 

7. R.C. Hontgomary and H.G. Hatch; 'Application of Differential 
Synthesis to Design of Multiaxis Stability Augmentation 
System', Journal of Aircraft, Vol. 6, Aug. 1969. 

8. J.P. Campbell and M.O. Mckinney; 'Summary of Methods for 
Calculating Dynamic Lateral Stability and Response and for 
Estimating Lateral Stability Derivatives' , BACA Report 
1098/52. 

9. G.¥. Hall: 'A. Flight Test Investigation of Direct Side Force 
Control', Air Force Flight Dynamics Lab,, TR-71-106 Sept. 

10. E. Yore: 'Optimal Decoupling Control', Proc. 1968, JACC, 

IEEE, Bew York, pp. 327-336. 



112 


Books 


11. y.A. yolovitcli; 'Linear Multivariable Systems', 

Sp r ing e r ~ Ve r lag . 

12. C.T. Chen; ' Intrcducxion to lineai’ System Theory', Eolt, 
Riiaeha.?rt and Winston. 


13. T.S. .Tortmann and K.L. liitz; 'An Introduction to linear 
Goiit.;-ol Systems' , Marcel Dckker 

14. '[illiaiii 1. Brogan. 'Modern Control Theory', Quentum, 

Hevj- York. 


15. 

16 . 

17. 

18. 


19. 


20 . 


21 . 


Jan Roskamt 'Blight Bynemics of Rigid and Elastic Air- 
planes' (Part One and Two), The University of Kansas, 
Lawrence, Kansas. 

Sqn.Ldr. B. Lickinson: 'Aircraft Stability end Control for 
Pilots and Engineers' , Sir Isaac Pitman and Sons Ltd. 


S. Seckel; 'Stability and Control of Airplanes and 
Helicopters', A.cademic Press, Hew York and London. 

Anon.i 'Plying Qualities of Piloted lirpla.nes' , Military 
Specifications MIl-F'-8785B(i.3G) , Sept. 1964. 


P.R. G-antmacheri 'The Theory of Matrices', Vol. 1 , 
Chelsea, Hew York. 

System/360 Scientific Subroutine Package ( 360A-CM-03K) 
Version III Programmer's Manual, Form H20-0205-O, l9o8, 
IBM Corporation. 


H. Coldstein: 'Classical Mechanics' 


A-ddison-Uesley . 



A70545 


&jE 


,,,p ,'^yt^ 



